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PREFACE 


This  report  describes  work  performed  on  atmospheric-radiative- 
transfer  models  as  part  of  the  Army's  target  signature  analysis  program. 
In  the  development  of  a comprehensive  model  which  embodies  reflectance, 
illumination,  transmittance,  and  detector  submodels  we  have  considered 
in  some  detail  an  illumination  model  for  poor  atmospheric  conditions. 

In  particular,  we  analyzed  mathematically  the  natural  radiation  field 
at  the  Earth's  surface  arising  from  the  moon  in  various  phases.  We  also 
examined  the  illumination  on  the  surface  as  a result  of  man-made  flares 
embedded  within  fogs . 

For  several  years  the  Environmental  Research  Institute  of  Michigan 
(ERIM)  has  been  developing  atmospheric-rad iative-transfer  models  pri- 
marily for  use  in  the  analysis  of  multispectrai  data  acquired  by 
satellites  and  aircraft.  The  present  study  extends  our  capability  to 
analyze  artificial  sources  and  extremely  dense  atmospheres. 

The  research  described  in  this  report  was  performed  under  Contract 
DAAA21-74-C-0331  and  covers  the  period  from  18  March  1974  through 
10  April  1975.  Mr.  Matthew  Nowak  served  as  Technical  Monitor.  The  pro- 
gram was  directed  by  R.  R.  Legault,  Vice  President  of  ERIM  and  R.  E. 
Turner,  Principal  Investigator.  The  ERIM  number  of  this  report  is 
108300-4-F. 
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D.  Dickerson  for  her  secretarial  assistance  in  the  preparation  of  this 
report . 
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1 

SUMMARY 


Over  the  last  seven  yearn  considerable  effort  has  been  spent  on 
atmospheric-radiative-transfer  modeling  at  the  Environmental  Research 
Institute  of  Michigan  (ERIM) . Much  of  this  effort  has  been  associated 
with  remote  sensing  applications  using  aircraft  and  satellite  multi- 
spectral  data.  This  report  is  a continuation  of  atmospheric-radiative- 
transfer  modeling  techniques  with  particular  emphasis  on  low  visibility 
atmospheres  and  artificial  sources. 

The  impetus  for  this  study  arises  from  a need  by  the  Army  to  develop 
a comprehensive  model  for  the  target  signature  analysis  program  which 
can  be  used  to  determine  illumination  levels  on  military  scenes  as  a 
result  of  natural  and  artificial  sources.  In  this  investigation  we 
considered  in  detail  the  natural  illumination  on  the  Earth's  surface  as 
a Tesult  of  the  moon  in  various  phases  and,  in  particular  for  various 
atmospheric  states.  It  was  found  that  the  illumination  is  dependent  upon 
atmospheric  visibility  and  surface  reflectance  in  a predictable  but  some- 
what surprising  way. 

A major  pan.  of  this  study  involved  the  analysis  of  the  illumination 
levels  due  to  one  or  more  flares  (isotropic  point  sources)  embedded 
either  within  or  outside  of  dense  fogs.  It  was  found  that  spatial 
distribution  over  the  surface  was  quite  extended  for  flares  above  the 
fog  layer  but  more  concentrated  when  the  flare  was  inside  the  fog.  Con- 
tours are  plotted  which  show  the  illumination  levels  for  several  flares 
in  various  configurations. 

In  the  final  section  of  this  report  we  discuss  the  advances  made  as 
a result  of  this  investigation  and  we  recommend  additional  studies  to 
continue  the  current  modeling  effort  in  order  to  develop  the  most 
realistic  atmospheric  radiation  for  target  signature  analysis. 
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2 

INTRODUCTION 


In  order  to  model  the  real  world  of  illumination  it  is  necessary 
that  experiments  be  conducted  to  determine  typical  values  of  quantities 
which  one  ultimately  would  like  to  measure  or  know  in  order  to  quantify 
a realistic  environment.  It  is  oftentimes  not  practical  to  conduct  the 
possibly  enormous  number  of  measurements  which  would  be  necessary  to 
define  all  environmental  conditions  however.  Fo.  this  reason,  it  is 
desirable  that  one  develop  a model  of  the  realistic  environment  and  test 
the  model  calculations  with  experimental  data.  Such  a procedure  will 
then  determine  the  range  of  validity  of  the  model  and  enable  the  user 
to  gain  confidence  in  the  use  of  model  calculations  to  simulate  specific 
missions . 

Although  a number  of  measurements  of  the  natural  nighttime  illumina- 
tion of  the  Earth's  surface  have  been  made  over  the  years,  there  have  been 
no  attempts  to  model  realistic  atmospheres  and  the  nighttime  illumination 
through  these  atmospheres  starting  from  the  basic  principles  of  radiative- 
fansfer.  In  subsequent  sections  of  this  report  we  will  make  use  of  an 
atmospheric-radiative-transfer  model  to  calculate  the  natural  nighttime 
spectral  irradiance  on  the  Earth's  surface  as  a result  of  the  moon.  It 
should  be  noted  that  the  procedure  of  developing  a mathematical  model  is 
generally  preferable  to  the  method  of  using  a strictly  empirical  model  of 
the  radiation  field.  This  is  true  because  the  great  variation  in  environ- 
mental factors  can  preclude  practical  use  of  empirical  models.  Thus, 
making  use  of  an  empirical  model  of,  say,  illumination  vs.  visibility, 
would  lead  one  to  consider  only  that  factor  when,  in  reality  other  con- 
ditions such  as  fog  profile,  e.g.,  might  have  a more  important  effect  on 
the  overall  results.  Therefore,  a mathematical  model  based  upon  sound 
physical  principles  should  be  of  greater  utility  to  the  user  in  the  prac- 
tical implementation  of  the  model  to  simulate  reality. 
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In  a later  section  we  will  develop  a model  of  artificial  point 
sources  (flares)  embedded  in  various  atmospheres  of  varying  degrees  of 
turbidity.  The  model  developed  is  approximate  only  in  the  sense  that 
multiple  scattering  is  neglected.  Realistic  atmospheric  states  may  be 
considered  and  the  results  are  valid  for  wavelengths  in  the  visible  and 
near  infrared,  except  for  the  bands  of  gaseo.s  absorption. 

Finally,  we  model  several  point  sources  located  at  specified  posi- 
tions in  a horizontally  uniform  atmosphere  which  is  also  being  illumi- 
nated by  moonlight.  Contours  of  the  irradiance  are  presented  which  show 
the  relative  magnitude  of  artificial  and  natural  radiation  fields . With 
these  computational  procedures  one  can  simulate  a great  variety  of  situa- 
tions involving  natural  and  artificial  sources  of  illumination  under 
nighttime  conditions. 
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3 

ATMOSPHERIC  STRUCTURE 


In  this  chapter  we  shall  discuss  the  physical  state  of  the  atmosphere 
in  general  terms.  Thus,  we  shall  characterize  the  weather  in  terms  of  the 
static  conditions  of  the  gaseous  component  and  the  particulate  component, 
i.e.,  we  exclude  specific  temporal  changes  due  to  wind  and  precipitation. 
We  shall  define  the  complete  optical  state  of  the  atmosphere  for  all  con- 
ditions in  which  precipitation  does  not  occur. 

3.1  GASES 

Our  atmosphere  is  composed  of  the  major  permanent  gaseous  components, 
molecular  oxygen,  nitrogen,  and  argon.  These  gases  absorb  very  little 
radiation  in  the  visible  and  near-infrared  part  of  the  spectrum,  whereas 
the  variable  gaseous  components,  water  vapor,  carbon  dioxide,  ozone, 
sulfur  dioxide,  nitrogen  compounds,  methane,  and  other  trace  gases  do 
absorb  a significant  amount  of  radiation.  Of  these  components,  ozone, 
water  vapor,  and  carbon  dioxide  are  the  most  important  absorbers  of  radia- 
tion in  the  near-infrared.  If  other  gases  are  introduced  to  the  immediate 
environment  where  observations  are  being  made  then  some  quantitative  para- 
meterization of  their  composition  and  concentration  must  be  made  in  order 
to  account  for  the  attenuation  of  radiation.  It  is  beyond  the  scope  of 
this  report  to  consider  gaseous  absorption  in  detail.  More  detailed 
studies  of  atmospheric  gases  have  been  done  by  Goody  [1], 

3.2  AEROSOLS 

One  may  define  an  aerosol  as  a semi-permanent  suspension  of  solid  or 
liquid  particles  in  the  Earth's  atmosphere.  Typical  aerosol  distributions 
would  be  hazes,  clouds,  mists,  fogs,  smokes,  and  dusts.  It  is  the  aero- 
sol component  of  the  atmosphere  which  is  quite  variable.  Depending  upon 
location  and  time,  aerosols  can  have  many  compositions,  sizes,  shapes, 
and  densities.  In  condensation  processes,  the  predominant  composition 

^R.  M.  Goody,  Atmospheric  Radiation,  Oxford  University  Press, 

New  York,  1964. 
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is  water  (for  which  the  refractive  index  is  well  known) . Depending  upon 
the  history  of  the  distribution,  however,  various  contaminants  can  mix 
with  water  and  water  vapor  in  the  formation  of  droplets  to  produce  a com- 
position quite  different  from  that  of  water.  The  major  primary  source  of 
atmospheric  aerosols  on  a world-wide  basis  is  sea  t>pray  [2];  other  primary 
sources  are  wind-blown  dust  and  smoke  from  forest  fires.  While  the  com- 
position of  a world-wide  aerosol  distribution  can  be  estimated,  this  com- 
position can  nevertheless  vary  over  a wide  range  depending  upon  its 
proximity  to  various  local  sources . The  only  ways  in  which  the  composi- 
tion can  be  determined  is  by  direct  in  situ  sampling  of  the  air  in  a 
particular  location  or  by  radiometric  techniques.  Much  in  situ  sampling 
has  been  done  by  Volz  [3,41,  Grams  et  al.  [5],  and  Flanigan  and  Delong  [6]. 
As  a result  of  their  work  and  that  of  others , the  mean  world-wide  aerosol 
can  be  considered  to  have  a real  refractive  index  of  about  1.5  and  an 
imaginary  part  between  0.01  and  1.0.  Hence, 

m(A)  = m1  - im2  (1) 

where  m^  and  m2>  respectively,  art  the  real  and  imaginary  part  of  the 
refractive  index,  and  A is  the  wavelength.  For  realistic  atmospheric 
conditions , 

1.33  < mx  < 1.55  (2) 


[2] 

G,  M.  Hidy  and  J.  R.  Brock,  An  Assessment  of  the  Global  Sources 
of  Tropospheric  Aerosols,  Academic  Press,  New  York,  1971,  p.  1088. 

[3] 

JF.  E.  Volz,  Infrared  Absorption  by  Atmospheric  Aerosol  Substances, 
J.  Geophys.  Res.,  77,  No.  6,  1972,  p.  1017. 

[4] 

F.  E.  Volz,  Infrared  Optical  Constants  of  Ammonium  Sulfate,  Sahara 
Dust,  Volcanic  Pumice,  and  Flyash,  Appl.  Optics,  12,  No.  3,  1973,  p.  564. 

^G.  W.  Grams,  I.  H.  Blifford,  Jr.,  E.  G.  Schuster,  and  J.  J. 
DeLuisi,  Complex  Index  of  Refraction  of  Airborne  Fly  Ash  Determined  by 
Laser  Radar  and  Collection  of  Particles  at  13  km,  J.  of  Atmos.  Sciences, 
29,  1972,  p.  900. 

^D.  F.  Flanigan  and  H.  P.  DeLong,  Spectral  Absorption  Characteris- 
tics of  the  Major  Components  of  Dust  Clouds,  Appl.  Optics,  10,  No.  1, 

1971,  p.  51. 
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and 


(3) 


0 < m2  < 1.0 

which  are  roughly  independent  of  wavelength. 

Particles  also  come  in  various  sizes.  Generally,  however,  they  can 
be  put  into  three  categories:  the  Aitken  nuclei  with  radii  between  10”  ^ 

and  10"5  cm,  the  large  particles  with  radii  between  10'"'’  and  10“^  cm,  and 
the  so-called  giant  particles  with  radii  greater  than  10“^  cm.  In  most 
hazes,  the  optically  active  region  is  made  up  of  those  particles  in  the 
large  or  giant  categories.  Junge  [7]  showed  that  most  aerosol  distribu- 
tions follow  the  simple  power  law: 

n(z,  r)  ■ C(z)r“v  (4) 

where 

n(z,  r)  is  the  particle  number  density  for  radius,  r, 
at  altitude,  z 

v is  the  exponent  of  the  power  law 

Experimentally,  v has  been  found  to  vary  from  2 to  5 for  various  tropo- 
spheric aerosol  distributions.  Also,  other  investigators  have  found  a 
better  fit  to  the  particulate  data  by  using  the  modified  gamma  distribu- 
tion, 

n(z,  r)  - ara  exp  (-brY);  0 <.  r < <*>  (5) 

where  a,  a,  b,  and  y are  parameters  which  describe  the  distribution. 

For  the  liquid  aerosols  it  can  be  assumed  that  the  particles  are 
spherical  or  nearly  spherical  in  shape.  For  solid  particles,  however, 
the  shape  may  assume  any  form.  A number  of  investigators  [8,9,10]  have 

™C.  E.  Junge,  Journal  Meteorol.,  12,  1955,  p.  13. 

^A.  C.  Holland  and  G.  Gagne,  The  Scattering  of  Polarized  Light  by 
Polydisperse  Systems  of  Irregular  Particles,  Appl.  Optics,  9,  No.  5,  1970, 
p.  1113. 

Wr.  S.  Powell,  R.  R.  Circle,  D.  C.  Vogel,  P.  D.  Woodson  III,  and 
B.  Bonn,  Optical  Scattering  from  Non-Spherical  Randomly  Aligned,  Poly- 
disoerse  Particles,  Planet.  Space  Sci.,  15,  1967,  p.  1641. 

^°^D.  M.  Ruthven  and  K.  F.  Loughlin,  Effects  of  Particle  Shape  and 
Size  Distribution  on  the  Transient  Solution  of  the  Diffusion  Equation, 
Nature  Phys.  Sci.,  230,  1971,  p.  69. 


12 


studied  the  influence  of  particle  shape  on  the  scattering  of  radiation. 

The  scattering  of  electromagnetic  radiation  by  odd-shaped  particles  has 
a different  pattern  from  that  produced  by  spherical  particles;  but  given 
a polydisperse  collection  of  odd-shaped  particles,  the  nature  of  the 
difference  is  not  clear.  Most  of  the  work  on  radiation  in  atmospheres 
concerns  spherical  aerosols;  we  shall  follow  suit  and  neglect  the  compli- 
cations of  particle  shape. 

3.3  OPTICAL  WEATHER 

Generally,  we  can  classify  almost  all  states  of  the  atmosphere  as 
being  either  hazy  or  foggy.  The  haze  may  be  almost  non-existent  (clear) 
or  very  heavy.  Fogs  begin  where  haze  or  mist  ends,  i.e.,  the  optical 
attenuation  is  usually  much  greater. 

3.3.1  HAZES 

Haze  can  be  defined  as  a state  of  the  atmosphere  when  the  visibility 
is  less  than  the  ultimate  visibility  of  336  km  for  an  aerosol-free  atmos- 
phere. For  the  casual  observer  the  usual  region  of  haze  is  that  condition 
for  which  the  visual  range  (at  sea  level)  in  a horizontal  path  lies 
between  23  km  and  1 km. 

Deirmendjian  [11]  has  used  the  modified  gamma  distribution  (Equation 
(5))  to  characterize  different  hazes.  The  following  is  a list  of  the 
hazes  Deirmendjian  considered  and  the  relevant  parameters. 


Haze  Type 

a 

a_ 

1 

b 

M 

5.3333  x 104 

1 

1/2 

8.9443 

L 

4.9757  x 106 

2 

1/2 

15.118b 

H 

4.0000  x 105 

2 

1 

20.0000 

Haze  M,  used  to  describe  a marine  or  coastal  haze,  has  a peak  in 
the  distribution  at  r - 0.05  pm.  The  haze  L represents  a continental 
distribution  and  has  a peak  at  r * 0.07  pm.  The  haze  H model  can  be  used 

^^D.  Deirmendjian,  Electromagnetic  Scattering  on  Spherical  Poly- 
dispersions, Elsevier,  New  York,  1969,  p.  78. 
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to  represent  a stratospheric  aerosol  or  dust  layers;  it  has  a peak  at 
0.10  um.  A graph  of  the  three  hazes  is  illustrated  in  Figure  1.  We  shall 
make  use  of  these  l.azes  in  the  detailed  analysis  of  atmospheric  radiation. 

Thus  far  we  have  considered  the  composition,  sizes,  and  shapes  of 
aerosol  particles,  but  in  order  to  define  the  atmospheric  state  we  must 
also  knew  the  concentration  of  particles.  Wiegand  [12]  was  th«  first  to 
measure  the  vertical  profile  of  an  aerosol  number  concentration  of  con- 
densation nuclei.  As  a result  of  his  measurements  and  many  others  cited 
in  Ivlev  [13] , it  was  found  that  the  concentration  of  condensation  nuclei 
obeys  an  exponential  law  with  altitude.  It  was  also  determined  that  a 
zone  of  increased  concentration  of  large  particles  exists  in  the  17  to 
23  km  altitude  range  (Junge  layer)  and  in  a probable  layer  under  the 
tropopause  at  9 to  10  km.  These  layers  are  relatively  stable  compared  to 
the  lower  part  of  the  troposphere.  As  an  example  of  an  aerosol  model 
atmosphere,  Zuyev  [14]  has  constructed  the  following  altitude-dependent, 
number  density-distribution-function  lor  aerosols: 


/N(0)e“bz;  z 1 5 km 
0.03;  5 km  <.  z £ 15  km 

10.03e0,06z;  15  km  £ z £ 20  km 

0.01e"0,09z;  z > 20  km 


(6) 


where  N(z)  and  N(0),  respectively,  are  the  number  densities  at  altitude  z 
and  the  Earth's  surface  (altitude  U) . 


[121 

C.  Wiegand,  Chemical  Composition  and  Radioactivity  of  the  Atmos- 
phere, C.  Junge,  World  Publishing  Co.,  Cleveland,  1963. 

[131 

L.  S.  Ivlev,  Aerosol  Model  of  the  Atmosphere,  Prob.  Fiz.  Atmos. 
Leningrad,  No.  7,  1967,  pp.  125-160. 

[14] 

V.  Ye.  Zuyev,  Propagation  of  Visible  and  Infrared  Waves  in  the 
Atmosphere,  NASA  TTF-707,  NTIS,  U.  S.  Department  of  Commerce,  Springfield, 
Va.,  1970. 
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PARTICLES  PPR  UNIT  VOLUME  n(r) 


FIGURE  1.  HAZE- TYPE  DISTRIBUTION  FUNCTIONS  USED 
Units  depend  on  the  particular  model. 


In  conclusion,  we  can  say  that  the  most  significant  aspect  of  aero- 
sol particles  in  the  atmosphere  is  their  high  degree  of  variability  — in 
composition,  size  distribution,  and  especially  in  number  density.  All  the 
models  in  the  current  literature  deal  with  a highly  approximate  average- 
condition  from  which  large  deviations  can  occur  in  real  situations.  The 
many  details  of  aerosol  science  will  not  be  considered  in  this  report. 

For  a more  complete  study  of  the  physics  and  chemistry  of  aerosols,  one 
should  consult  Mason  [15],  Fuchs  [16],  Davies  [17],  or  Green  and  Lane  [18]. 

3.3.2  FOGS 

Fo*;,-  represent  those  conditions  in  which  the  visual  range  ib  less 
than  1 km.  Thus,  in  the  visible  part  of  the  spectrum  there  can  be  severe 
attenuation  of  radiation.  In  contrast  to  hazes  the  relative  humidity  in 
a fog  is  very  high  and  condensation  of  water  vapor  onto  small  nuclei  has 
taken  place.  Fogs  and  clouds  are  quite  similar  in  terms  of  their  optical 
properties  but  the  conditions  under  which  they  form  may  be  quite  different. 
To  understand  the  conditions  for  fog  formation  let  us  refer  to  figure  2. 

The  curve  in  Figure  2 represents  the  condition  in  which  the  air  is  just 
saturated  with  water  vapor.  At  a given  temperature  it  will  hold  the 
amount  indicated  by  the  curve.  This  is  called  the  dew  point.  If  the 
temperature  falls  then  the  air  is  supersaturated,  condensation  takes  place 
and  a fog  or  cloud  forms . As  a consequence , the  dew  point  decreases 
because  water  vapor  has  condensed . 

Fogs  are  generally  classified  according  to  the  conditions  under  which 
they  form.  There  are  three  general  categories:  air  mass  fogs,  frontal 

fogs,  and  ice  fogs. 

J.  Mason,  The  Physics  of  Clouds,  Clarendon  Press,  Oxford,  1971. 

A.  Fuchs,  The  Mechanics  of  Aerosols,  Pergaramon,  Elmsford,  N.Y., 

1964. 

117>C.  N.  Davies,  Aerosol  Science,  Academic  Press,  New  York,  1964. 
f 181 

1 JH.  L.  Green  and  W.  R.  Lane,  Particulate  Clouds,  Dusts,  Smokes, 
and  Mists,  2nd  Edition,  D.  Van  Nostrand,  Inc.,  Princeton,  1964. 
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TEMPERATURE  IN  DECREES  FAHRENHEIT 


FIGURE  2.  THE  SATURATION  CURVE 


Air  mass  fogs  can  be  subdivided  into  four  general  types:  1)  advec- 

tion,  2)  radiation  fogs,  3)  advection-radiation  fogs,  and  4)  upslope  fogs. 
Advection  fogs  form  when  warmer  air  blows  over  a colder  land  or  water 
surface.  Thus,  the  warmer  air  gives  off  heat  to  the  colder  surface, 
thereby  cooling  the  air  to  the  dew  point.  The  water  vapor  in  the  air 
condenses  on  the  nuclei  of  salt,  dust,  or  smoke  even  before  the  relative 
humidity  reaches  100  per  cent.  The  air  does  not  tend  to  rise  because  it 
is  being  cooled  and  it  then  becomes  dense,  heavy,  and  stable.  Examples 
of  advection  fogs  are  the  following: 

1.  Land-  and  Sea-Breeze  Fog:  Characterized  by  warm,  moist,  land  air 

passing  over  cold  water.  Typical  of  New  England  coast. 

2.  Sea  Fog:  This  is  similar  to  the  land-  and  sea-breeze  fog  except 
that  it  arises  from  the  cooling  of  sea  air  over  a cold  ocean  current  or 
cold  coastal  current.  Typical  of  California  coast  in  summer. 
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3.  Tropical-Air  Foy:  This  is  caused  by  the  gradual  cooling  of 

tropical  air  as  it  moves  poleward  over  the  ocean  or  over  land.  It  is 
probably  the  most  coonon  type  of  fog  over  the  open  sea.  Typical  of  the 
southeastern  U.S.  and  east  coast  of  the  U.S. 

4.  Steam  Fog:  Steam  fog,  characterized  by  tufts  of  whirling  smoke 

coming  cut  of  the  water  is  caused  by  the  passage  of  cold  air  over  warmer 
water.  In  polar  latitudes  it  is  called  arctic  sea  smoke. 

Radiation  fogs  are  different  from  advection  fogs  in  that  they  form 
over  land  when  the  sky  is  clear,  the  wind  is  light,  and  the  relative 
humidity  is  high  near  sunset.  Since  the  relative  humidity  is  high  a 
slight  cooling  will  result  in  a lowering  of  the  temperature  to  the  dew 
point.  We  have  two  kinds  of  radiation  fogs. 

1.  Ground  Fog:  Tills  is  the  simplest  type,  involving  light  winds, 

clear  skies,  and  a surface  temperature  inversion.  Typical  of  the 
Appalachian  valleys. 

2.  High  Inversion  Fog:  A winter  time,  land  phenomenon  resulting 

from  more  than  a single  night  of  cooling.  This  is  typical  of  the 
tropical  west  coasts  of  continents  as  in  California. 

The  third  type  of  alt  mass  fog  is  the  advection-radlation  fog.  It 
is  caused  by  nighttime  radiational  cooling  of  air  that  has  moved  inland 
from  the  sea  during  the  day.  Typical  of  the  California  coast  or  the 
Great  Lakes  during  Fall. 

The  fourth  type  of  air  mass  fog  is  the  upslope  or  expansion  fog. 

This  Involves  the  cooling  of  air  in  its  dew  point  by  adiabatic  expansion 
as  the  air  moves  to  higher  elevations.  It  is  typical  of  the  Great  Plains 
in  the  U.S. 

There  are  three  kinds  of  frontal  fogs,  prefrontal,  post-frontal, 
and  front-passage  fogs. 
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1.  Prefrontal  Fog  (warn  front):  This  is  characterized  by  rain 

falling  into  stable  cooler  air,  evaporating  and  raising  the  dew  point 
of  the  air  to  the  condensation  point.  It  is  typical  of  the  Mid-  and 
North-Atlantic  coast  states  in  Winter;  also  the  Appalachian  valleys. 

2.  Post-Frontal  (cold  front):  This  fog  is  similar  to  the  pre- 

frontal fog  except  that  the  associated  precipitation  band  .s  much  more 
restricted  in  area.  It  is  typical  of  the  Midwestern  U.S.  during  polar 
air  mass  outbreaks . 

3.  Front-Passage  Fog:  This  is  a short  lived  fog  that  can  form  in 

various  ways  with  the  passage  of  a frontal  zone. 

Finally,  we  have  the  last  category  of  fogs,  the  ice  fog.  It  is 
caused  primarily  by  man-made  sources  of  moisture  (engines  runnings)  at 
temperatures  lower  than  about  -30°C.  Typical  of  the  polar  regions. 

Having  considered  the  general  types  of  fogs  and  their  formation  we 
shall  now  investigate  the  size  distribution  of  the  droplets.  When  drop- 
lets begin  to  grow  in  size  at  the  saturation  point  the  growth  proceeds 
at  a regular  rate.  Neiburger  and  Chien  [19]  investigated  this  growth 
process  and  found  that  the  size  distribution  of  the  particles  changes  as 
the  fog  ages  and  that  the  distribution  becomes  bimodal.  This  is  illus- 
trated in  Figure  3.  As  time  goes  on  and  saturation  and  then  super- 
saturation occurs  the  transition  from  haze  to  fog  is  quite  rapid.  As 
condensation  takes  place  the  particles  coalesce  and  the  fog  siz«  distri- 
bution develops  a secondary  peak  in  the  3 co  6 wm  range. 

The  shape  of  the  particles  are  assumed  to  be  spherical  since  in  a fog 
we  have  water  droplets.  Their  composition,  however,  may  vary  since  a fog 
may  be  contaminated  by  dust,  smoke,  or  other  matter  in  the  immediate 
environment . 


[19] 

M.  Neiburger  and  C.  W.  Chien,  Physics  of  Precipitation,  Geophys. 
Monogr.  No.  5,  Washington,  D.C.,  1960,  pp.  191-210. 
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CONCENTRATION  (CM'3)  PEN  MICRON  RAOIUS 


FIGURE  3.  AEROSOL  DISTRIBUTIONS  AS  A FUNCTION  OF  I, APSE 
TIME  AFTER  THE  INITIATION  OF  THE  CONDENSATION 
PROCESS.  The  times,  in  minutes,  are  indicated 
at  the  base  of  each  curve  after  Neiburger  and 
Chien  (1960). 
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4 

ATTENUATION  OF  RADIATION 


In  this  section  we  shall  investigate  the  interaction  of  visible  and 
infrared  radiation  with  a scattering  and  absorbing  medium  such  as  a haze 
or  fog.  Before  a radiative-transfer  model  is  developed  we  need  to  know 
the  basic  input  parameters  for  the  environment  under  consideration. 

4.1  TRANSMITTANCE 

The  transmittance  of  a beam  of  radiation  propagating  through  the 
atmosphere  lies  somewhere  between  zero  and  unity.  For  zero  transmittance 
no  direct  radiation  passes  through  the  medium  and  for  a value  of  one  all 
the  radiation  passes  through  undiminished  in  intensity.  We  shall  neglect 
gaseous  absorption  in  this  report  because  we  are  primarily  interested  in 
the  scattering  of  radiation  by  gases  and  aerosols . 

If  an  infinitesimal,  narrow  beam  of  monochromatic  radiation  propagates 
through  a medium  the  intensity  at  a distance  x from  the  source  is 

IA(x)  = Ix(o)Ta(jO  (7) 

where  I^(o)  is  the  source  intensity  at  wavelength  X and  T^ (x)  is  the 
spectral  transmittance  at  distance  x.  If  the  loss  of  radiation  from  the 
beam  is  total,  i.e.,  complete  absorption  and/or  scattering  out  of  the  beam 
then  the  fractional  change  in  intensity  is  given  by  the  following: 

AIa(x) 

1^00"  ' KJ(X)  <8) 

where  k^(x)  is  the  probability  par  unit  length  for  a net  loss  of  radiation 
from  the  beam  and  Ax  is  the  distance  over  which  the  fractional  change 
occurs.  As  the  length  Ax  approaches  zero  we  have 

dl,(x) 

~dT='^(x)  VX)  (9) 
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or 


K^(x')dx' 


Ix(x)  - Ix(o)  e 


(10) 


Thus,  we  see  that  the  spectral  transmittance  T, (x)  is  given  by 


Kx(x')dx' 


Tx(x)  = e 


(ID 


Hence,  one  of  the  problems  is  to  determine  the  probability  per  unit  length 
<x(x) . 

4.2  ATTENUATION  COEFFICIENTS 

The  probability  per  unit  length  <x(x)  is  also  called  the  volume 
extinction  coefficient  of  the  medium.  In  this  section  we  shall  see  how 
this  quantity  is  determined  for  various  atmospheric  states. 

4.2.1  VISIBILITY 

The  spectral  transmittance  of  a region  of  the  atmosphere  can  be  deter- 
mined oy  measuring  the  intensity  Ix(x)  of  a laser  beam  at  distance  x from 
the  laser.  There  are  times  however,  when  it  is  not  convenient  to  use 
lasers  in  field  measurements.  We  would  then  like  to  have  another  method 
of  estimating  the  extinction  coefficient  k^(x). 

Equation  (10)  describes  the  intensity  of  a self-luminous  object  at 
distance  x.  If,  in  addition  there  is  background  radiation  being  received 
by  a detector  at  distance  x then  we  have 


!x(x) 


Ix(o)  T(x)  + 


Bx(x')tcx(x')T(x,)dx’ 


o 


(12) 
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where  the  integral  represents  the  sum  of  all  background  radiation  along  a 
given  path.  These  components  are  depicted  in  Figure  4.  If  the  background 
radiation  B^x’)  is  independent  of  distance  as  it  usually  is  in  homo- 
geneous atmospheres  then  Equation  (12)  becomes 

Vx>  " V °)T(x)  + Ba[1-T(x)]  (13) 

Expression  (13)  is  quite  general  and  applies  to  a variable  path.  If  the 
luminous  object  is  moved  to  a great  distance  then  the  intensity  is  just 
about  due  to  the  atmosphere,  i.e., 

Vx)  * BA  as  x - 00  (14) 

If  the  background  is  defined  as  the  atmosphere  then  the  contrast  can  be 
defined  as 


c - ia(o)  ~ Vx)  _ (yp)  - y [i-t(x)  3 
Vx)  Ia(o)T(x)  + BAfl-T(x)J 

Now,  if  an  intrinsic  contrast  Cq  is  defined  as 


C 

o 


Vo)  - BA 


(15) 


(16) 


then 


C = C T 
o c 


(17) 


where  the  contrast  transmittance  T is 

c 


T 

c 


B.  l-T(x) 


(18) 


Thus,  by  measuring  the  contrast  the  distance 
case  of  a black  object  measured  against  the 
to  Equation  (13) 


x cai  be  determined.  For  the 
atmosphere  we  have  according 


IA(x)  = Ba[1-T(x)] 


(19) 


23 


TARGET 


DETECTOR 


FIGURE  4.  DETECTION  OF  DIRECT  AND  BACKGROUND  RADIATION 


Thus,  the  contrast  Is 


C 


Vx) 


(20) 


or 


k(x* Jdx 


C = T(x) 


(21) 


Koschmieder  [20]  defined  the  limiting  contrast  for  the  eye  as  0.02.  In 
that  case  the  distance  at  which  the  contrast  is  0.02  for  a constant  atmos- 
phere is 

V = — -1"  (22) 

K 


where  V is  the  so-called  visual  range  and  k is  the  volume  extinction 
coefficient  at  a wavelength  of  0.55  um.  Figure  5 illustrates  this  varia- 
tion of  visual  range  with  k . Hie  ultimate  visual  range  is  that  for  an 
aerosol  free  atmosphere.  This  value  is  336  km.  More  realistic  cases  are 
those  for  smaller  visual  ranges.  The  fog  regime  begins  rather  abruptly 
near  a visual  range  of  1 kilometer.  As  we  shall  see  in  a later  section  a 
mist  regime  may  also  be  defined  in-between  the  haze  and  fog. 

4.2.2  HAZES 

Knowing  the  gaseous  composition  of  the  atmosphere  one  can  calculate 
the  Rayleigh  extinction  coefficient  at  any  wavelength.  This  was  done  by 
one  of  the  authors  (Turner)  for  any  visible  and  near-infrared  wavelength 
and  as  a function  of  altitude.  A correction  was  also  made  for  variable 
terrain  altitudes.  The  Rayleigh  extinction  coefficient  is 

kr(z)  * oR  Nr(z)  (23) 

Koschmieder,  Theorie  der  horizontalen  Sichtweite,  Beitr. 

Phys.  freien  Atm.  12,  1924,  pp.  33-53. 
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EXTINCTION  COEFFICIENT  [*(km  J)at  0.55  4m] 


FIGURE  5.  VARIATION  OF  VISUAL  RANGE  WITH  EXTINCTION  COEFFICIENT 


where  o is  the  Rayleigh  cross  section  and  N (z)  is  the  molecular  number 
density.  The  cross  section  has  the  well-known  l/X4  wavelength  dependence 
For  a haze  there  is  more  scattering  but  since  the  scattering  centers  are 
large  compared  to  the  wavelength,  many  orders  of  scattering  contribute  to 
the  total  and  there  is  little  dependence  on  wavelength.  For  hazes  the 
spectral  dependence  is  approximately  X~l*3.  Figure  6 illustrates  the 
spectral  dependence  of  the  extinction  coefficient  for  hazes  and  clouds. 
These  coefficients  were  calculated  by  Deirmendjian  [11]  for  water  spheres 

The  dependence  on  altitude  is  more  complicated.  Elterman  [21] 
collected  data  on  various  experiments  to  determine  the  extinction  coef- 
ficients for  hazy  atmospheres.  The  results  are  presented  in  Figure  7. 
Here  we  see  a gradual  decrease  with  altitude  for  most  of  the  experiments. 
Elterman  [22]  has  also  modeled  hazes  in  terms  of  altitude,  visual  range, 
and  wavelength  making  use  of  experimental  data  at  a known  visual  range. 
The  resulting  compilation  is  what  we  shall  use  in  the  calculation  of  the 
natural  radiation  field. 

Hence,  we  can  say  that  the  model  described  above  characterizes  a 
general  "standard"  haze  but  there  can  be  large  variations  from  the  m,  del 
in  individual  cases.  In  a specific  situation  one  should  measure  either 
the  spectral  transmittance  or  the  visual  ranges  in  order  to  define  the 
state  of  the  atmosphere. 

4.2.3  FOGS 

Clouds  and  fogs  are  made  up  of  particles  which  are  much  larger  than 
those  for  hazes.  As  a result,  the  scattering  is  usually  less  selective, 
or  dependent  upon  wavelength.  Arnulf  et  al.  [23]  presented  data  on  fogs 

Elterman,  U.V.,  Visible,  and  IR  Attenuation  for  Altitudes  to 
50  km,  1968,  AFCRL  Report  No.  68-0153,  Office  of  Aerospace  Research. 

Elterman,  Vertical-Attenuation  Model  with  Eight  Surface 
Meteorological  Ranges  2 to  13  Kilometers,  AFCRL  Report  No.  70-0200,  1970, 
Office  of  Aerospace  Research. 

^^A.  Arnulf,  J.  Bricard,  E.  Cure,  and  C.  Veret,  Transmission  by 
Haze  and  Fog  in  the  Spectral  Region  0.35  to  10  Microns,  J.  of  the  Opt. 
Soc.  of  Amer.,  Vol.  47,  No.  6,  1957,  pp.  491-498. 
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ATTENUATION/SCATTERING  COEFFICIENT  IknT' 


WAVELENGTH  |jim) 


FIGURE  6.  SPECTRAL  ATTENUATION  AND  SCATTERING  COEFFICIENTS 
FOR  VARIOUS  HAZE  (L,M,H)  AND  CLOUD  (C.,C-,C  ) 
MODELS  (ADAPTED  FROM  DEIRMENDJIAN) . 1 6 
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AEROSOL  ATTENUATION  COEFFICIENT  (km”') 


ALTITUDE  (km) 


FIGURE  7.  AEROSOL  ATTENUATION  COEFFICIENTS  VS  ALTITUDE  AT  X = 0.55U 


Comparison  of  results: 

XXX X solar  aureole,  2 balloon  flights,  Newkirk  and  Eddy  (1964); 
a a a a solar  radiation  measured  from  aircraft,  mean  of  8 flights, 
Penndorf (1954); 

searchlight  probing,  mean  of  105  profiles,  Elterman 
(1966a)  and  (1966b}; 

• • • • balloon  integrating  nephelometer,  mean  of  14  flights,  Crosby 

and  Koerber (1962); 

• • • • solar  radiation  measured  from  balloons,  mean  of  3 flights, 

Kondratiev  et  al.  (1967); 

a a a a spacecraft  horizon  photography,  analysis  of  4 frames, 
Rosenberg  (1966),  Feoktistov  (1965); 
o o o o searchlight  probing,  mean  for  5 nights,  Rosenberg  (1966) 
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of  various  types.  The  extinction  coefficients  (density  per  kilometer) 
for  four  fogs  is  illustrated  in  Figures  8,  9,  10,  and  11.  It  should  be 
noted  that  selective  fogs  are  similar  to  hazes  in  that  there  is  a definite 
spectral  dependence  whereas  for  the  stable  fogs  there  is  essentially  no 
spectral  dependence  in  the  visible  and  near  infrared  regions,  in  fact, 
the  transmission  diminishes  in  the  near  infrared  in  some  cases.  Only 
oeyond  3.7  pm  does  the  extinction  coefficient  begin  to  decrease.  By 
measuring  droplets  in  a fog  and  comparing  these  measurements  with  calcu- 
lations for  various  fog  droplet  size  distributions  Eldridge  [24]  deter- 
mined four  basic  fog  extinction  coefficients.  These  are  depicted  in 
Figure  12.  Hence,  we  can  summarize  the  spectral  dependence  of  all  the 
coefficients  in  Figure  13.  These  results  are  taken  from  Deirmendjian  [25]. 
The  greatest  change  is  for  a pure,  Rayleigh  atmosphere,  and  as  the  medium 
develops  a haze  the  spectral  dependence  is  much  less  pronounced.  For  a 
fog  there  is  little  dependence  on  wavelength  in  the  visible  and  near 
infrared. 

In  order  to  calculate  the  volume  extinction  coefficient  we  must 
realize  that  the  total  coefficient  is  made  up  of  two  parts,  that  due  to 
the  gaseous  component  are  that  due  to  the  particles;  thus 

k(X,  r)  = <R(X,r)  + «A(X,m,r)  (24) 

where  <R(X,r)  is  the  Rayleigh  coefficient  at  position  in  space  r and  for 
wavelength  X and  <A(X,m,r)  is  the  aerosol  coefficient  at  the  same  point 
and  same  wavelength  but  for  a composition  denoted  by  the  refractive  index 
m which  is  wavelength  dependent . As  we  have  seen  the  Rayleigh  component 
is  easy  to  calculate  since  it  only  depends  on  the  gaseous  density  of  the 
atmosphere  and  the  wavelength.  The  aerosol  extinction  coefficient  is 

R.  G.  Eldridge,  Haze  and  Fog  Aerosol  Distributions,  J.  of  Atmos. 
Sciences,  1966,  pp.  605-613. 

f25l 

D.  Deirmendjian,  Quart.  Journal  Royal  Meteorol.  Soc.  86,  371, 

(1960) . 
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WAVELENGTH  (micrometers) 

FIGURE  8.  ATMOSPHERIC  TRANSMITTANCE  SPECTRA  FOR  SELECTIVE  FOGS 
AS  MEASURED  BY  ARNULF  ET  AL. 
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FIGURE  11.  ATMOSPHERIC  TRANSMITTANCE  SPECTRA  FOR  STABLE  FOGS, 
TYPE  2,  AS  MEASURED  BY  ARNULF  ET  AL. 
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ATTENUAT'ON  ( KM’1) 


FIGURE  12.  MEASURED  AND  COMPUTED  FOG  SPECTRAL  ATTENUATIONS 
ILLUSTRATING  GOODNESS  OF  FIT 


FIGURE  13.  TYPICAL  VALUES  FOUND  FOR  ATMOSPHERIC  EXTINCTION 
COEFFICIENT  VS  WAVELENGTH.  Indicated  values  for 
curves  B and  E are  P(180°)/4ir*Bgca. 
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however,  more  difficult  to  calculate.  If  uniform  spherical  particles  are 
assumed  with  a known  refractive  index  then  we  can  calculate  the  cross 
sections.  The  scattering  of  electromagnetic  radiation  by  homogeneous 
dielectric  spheres  is  called  MIE  scattering.  In  principle,  determining 
the  scattering  and  absorption  cross-sections  and  the  scattering  phase 
function  is  straightforward.  First,  the  wave  equation  is  solved  inside 
and  outside  the  sphere;  then  the  solutions  are  match  at  the  boundary  to 
determine  the  constants . The  intensity  of  the  radiation  scattered  into 
an  angle  x is  given  by 


I 


(25) 


where  R is  the  distance  from  the  sphere,  and  I is  the  incident  intensity. 

A • ^ 

The  scattering  amplitudes  | j z and  | S ^ ! are  given  by 


S 


S 


1 


2 


v 21+1  r . . , 

/.  i(i  + i)  ^Yi  + biTi^ 

i»i 

£(£  + D a*T*  + b£  £ 


(26) 


In  Equation  (26)  a£  and  b£,  which  are  the  Mie  coefficients,  are  given  by 


<J^(mx)<J>£(x)  - w}»£(mx)^(x) 
^(mx)c£(x)  - uu^raxKj^x) 


(27) 


nu|^(mx)»l>£(x)  - <Jj£(mx)t(>‘(x) 
mii>^(mx)c£(x)  - <|>£(mx)i^(x) 


(28) 
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where  m and  x,  respectively,  are  the  complex  refractive  index  and  size 
parameters  — that  is,  m ■ m^  - in^  and  x ■ 2^r/A  when  r is  the  particle 
radius.  The  functions  1 and  are  the  Riccati-Bessel  functions,  and  a 
prime  indicates  differentiation  with  respect  to  the  argument.  The  it  and 
functions  are  given  by 


and 


it£(cos  x) 


dPg,(cos  X) 
d cos  x 


(29) 


dTT£(cos  X) 

t£(cos  x)  " COS  X tt£(cos  x)  - sin  X d CQ8  — (30) 


where  P^(cos  x)  is  a Legendre  polynomial.  The  scattering  cross-section 
is 


Os(A,r,m) 


rtr  Q (x,m) 
s 


+ ib.n 


(31) 


where  Q (x,m)  is  called  the  scattering  efficiency  factor.  The  total 
8 

(scattering  plus  absorption)  cross-section  is 

9 X2  ® 

ot(A,r,m)  - nrzQt(x,m)  - — l (24  + l)Re(a£  + b£)  (32) 

where  Qt(x,m)  is  the  total  efficiency  factor.  Likewise,  the  absorption 
cross-section  is  then 

o (A,r,m)  - a (A,r,m)  - a (A,r,m)  (33) 

a L 8 

with  a corresponding  efficiency  factor  Q (x,m) . As  in  the  Rayleigh  case, 

& 

a scattering  phase  function  can  be  defined: 

— T Us,!2  + Is2I2J 

2ttx  Q (x,m) 


P(cos  x) 


(34) 
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Calculating  all  the  functions  given  above  for  various  refractive 
indices  and  size  parameters  is  quite  involved.  Nevertheless,  computer 
programs  have  been  written  which  allow  the  performance  of  this  analysis. 
Having  obtained  a program  written  by  J.  V.  Dave,  we  used  it  to  calculate 
the  cross-sections  and  phase  functions.  For  example,  we  have  calculated 
the  scattering  efficiency  factor  for  homogeneous  spheres  of  refractive 
indices  m = 1.29,  1.29  - 0.0465i,  and  1.28  - 1.37i.  The  results  are 
shown  in  Figure  14.  It  should  be  noted  that  the  efficiency  factor  is 
greatest  for  the  real  index  (that  is,  when  there  is  no  absorption)  and 
decreases  with  increasing  imaginary  index.  The  absorption  efficiency 
factor  has  also  been  calculated  for  m * 1.28  - 0.04651  and  1.28  - 1.37i. 
Here  the  efficiency  factor  is  large  for  a hi^h  imaginary  index  and  small 
particles,  as  illustrated  in  Figure  15.  Finally,  the  total  efficienc" 
factor  is  shown  in  Figure  16  for  the  same  set  of  refractive  indices. 

Thus,  it  can  be  seen  that  efficiency  factors  vary  strongly  for  different 
refractive  indices  and  size  parameters. 

Knowing  the  cross-sections,  one  can  then  calculate  the  scattering, 
absorption,  and  extinction  coefficients  by  multiplying  the  cross-sections 
by  the  particle  number  density.  For  real  atmospheric  conditions  charac- 
terized by  haze,  fogs,  and  dusts,  there  is  a distribution  of  particle 
sizes.  (One  distribution,  characterized  in  Section  3.2,  is  the  modified 
gamma  distribution.)  Thus,  for  a polydispersiou,  one  must  integrate  over 
particle  size  to  obtain  the  absorption,  scattering,  and  extinction 
coefficients: 


«A(*,  m,  z)  * [ A(X,  m,  r)n(z,  r)dr.  (35) 

0 


VA* 


z) 


O A<*.  m»  r)n<z»  r>dr 


(36) 
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SIZE  PARAMETER  (x  = 2irr/X) 


FIGURE  14.  SCATTERING  EFFICIENCY  FACTOR  FOR  HOMOGENEOUS 
SPHERES  OF  COMPLEX  REFRACTIVE  INDEX  m 
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FIGURE  15. 


ABSORPTION  EFFICIENCY  FACTOR  FOR  HOMOGENEOUS 
SPHERES  OF  COMPLEX  REFRACTIVE  INDEX  m 


FIGURE  16.  TOTAL  EFFICIENCY  FACTOR  FOR  HOMOGENEOUS 
SPHERES  OF  COMPLEX  REFRACTIVE  INDEX  m 
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<A(X,  tn,  z)  - I 0t  A(At  m,  r)n(z,  r)dr  (37) 

0 

where  a,  6,  and  < denote  absorption,  scattering,  and  extinction  coefficients 
A designates  aerosol 

n(z,  r)  is  the  aerosol-particle  number  density  at  altitude  z for 
particles  in  size  range  Ar  at  size  r 

The  number  density  is  normalized  as  follows: 


N(z)  = ( n(z,  r)dr 

i 

u 


(38) 


where  N(z)  is  just  the  total  number  density.  Likewise,  the  corresponding 
coefficients  for  molecular  scattering  can  be  found.  The  scattering  cross- 
section  is  well  known  and  the  absorption  cross-section  is  usually  taken 
to  be  that  for  ozone  in  the  visible  spectral  region.  Thus,  we  have  for 
the  complete  atmosphere: 

a(A,  m,  s,  z)  - aR(A,  z)  + aA(A,  m,  s,  z)  (39) 

B(A,  m,  s,  z)  - SR( A , z)  + 8a(A,  m,  s,  z)  (40) 

k(A,  m,  s,  z)  « <r(A,  z)  + <A(>,  m,  s,  z)  (41) 


where  we  have  shown  the  explicit  dependence  on  a complex  refractive  index 
m and  a particular  size  distribution  s. 

It  is  sometimes  useful  to  deal  with  an  average  cross-section  for 
aerosols.  Assuming  that  the  size  distribution  is  altitude-independent, 
we  can  write  Equations  (35),  (36),  and  (37)  as 

«A(*,  «,  z)  • N(z)“a>A(X’  (*2) 

Ba(A,  m,  z)  - N(z)os>A(A,  m)  (43) 
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where 


*A(*»  «.  z)  • N(z)ot^(X,  m) 


(44) 


°i,A(X’  n)  “ j °i,A(X*  “*  r)l^r)dr  (45) 

0 

In  which  the  function  4>(r)  is  normalized  to  1 and  the  index  i can  indicate 
a,  8,  or  t.  We  have  computed  many  values  of  the  average  cross-sections 
for  homogeneous  spheres  of  various  refractive  indices  and  size  distribu- 
tions. The  results  of  this  analysis  will  be  presented  later. 

The  angular  scattering  properties  of  the  spherical  particles  can  also 
be  calculated  using  Mie  theory.  We  will  not  go  into  the  details  of  this 
analysis  here.  There  are  several  excellent  treatises  on  this  subject, 
van  de  Hulst  [26],  Kerker  [27],  and  Deirmendjian  [11].  A typical  phase 
function  or  differential  scattering  cross  section  is  depicted  in  Figure  17 
for  the  case  of  a polydisperse  collection  of  cloud  or  fog  particles.  As 
one  can  see  the  radiation  is  sharply  peaked  in  the  forward  direction.  It 
is  this  anisotropy  which  presents  great  difficulties  in  the  mathematical 
analysis  of  radiative-transfer. 

Since  fogs  are  composed  of  water  droplets  it  seems  natural  to  relate 
the  visual  range  to  the  total  water  content  of  a fog.  It  was  noted  by 
Trabert  [28]  in  1901  that  the  visual  range  is  related  to  water  content  w 
by  the  simple  relation 

V ■ Cr/w  (46) 

C.  van  de  Hulst,  Light  Scattering  by  Small  Particles,  John 
Wiley  and  Sons,  1957. 

f271 

1 M.  Kerker,  The  Scattering  of  Light  and  Other  Electromagnetic 
Radiation,  Academic  Press,  1969. 

^®^W.  Trabert,  Die  Extintion  des  Lichtes  in  einem  truben  Medium 
(Sichtweite  in  Wolken) , Meteor.  Z.,  18,  pp.  518-525,  1901. 


44 


5 


where  C is  a constant  and  r is  the  droplet  radius.  This  follows  from  the 
fact  that  the  extinction  coefficient  < is 


V IT 


•i4,  n/\  -\ 


and  the  visual  range  is 


Vmhl U 

K\ 


(47) 


(48) 


where  1 * 0.55  v»m.  Now  from  Mie  theory  we  saw  that  for  very  large 
c 

particles  Q(X ,r)  ■+•  2 and  hence  is  independent  of  wavelength.  Thus , 


y - 3^912_ 
2 it  Nr2 


(49) 


Now,  the  water  content  w is  given  by 
4 3 

w - y tt  r Np 


(50) 


where  p is  the  density  of  a water  droplet.  Thus,  we  have 
V - 2.608  r/w 


(51) 


Therefore,  if  the  size  of  the  droplets  remains  constant,  the  visual  range 

varies  inversely  as  the  water  content.  If  however,  the  number  of  parti- 

1/3 

cles,  N remains  constant  then  we  have  r « w and  hence 

V - C'  w"2/3  (52) 

Instead  of  making  these  approximations  we  will  develop  a semi- 
empirical  model  of  fogs.  Consider  the  fog  extinction  coefficient 

00 

k(z,X>  * | n(z,r)  o(z,r,X)dr  (53) 

o 
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but.  If  o(z, r,A)  ■ o(r,A)  then 

o(r,A)  - trr2  Q(r,A)  (54) 

Furthermore,  let  the  size  distribution  be  given  by 

n(z,r)  * N(z)ii»(r)  (55) 

where  N(z)  is  the  number  of  particles  per  unit  volume  and  0(r)  is  the 
distribution  function.  Hence, 

00 

k(z,A)  = N(z)  [ \p(r)  vr2  Q(r,A)dr 

4 

0 

- N(z)  F(A)  (56) 


Now,  the  water  content  w(z)  is 


00 

w(z)  « [ n(z,r)  j irr3p  dr 
o 


- N(z)C 


Therefore,  we  have  the  scaling  relation 


<(z,A) 


"k  k 

k(z,A  ) K ( Zp , A ) 


K (Zp , A*) 


„ HiL.  wLe*_ 

F( A*)  w<zp> 


* f(^)  g(z)  K(zp,A*) 


(57) 


(58) 


(59) 


(60) 
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where  r Is  some  altitude  at  which  the  visual  range  Is  measured.  Thus, 

P 

we  have  a formula  to  describe  the  spectral  and  altitude  dependence  of 

ft 

the  fog  extinction  coefficient.  In  Equation  (60)  X >*0.55.  We  finally 
have  the  relation 


3.912  f (X)  g (z) 

Kij<z*x> nT) — 3 — <6i> 

We  shall  use  for  f^(X)  the  spectral  dependence  of  each  fog  type  1 
given  In  Figure  12.  The  function  8j(z)  which  defines  the  fog  profile  Is 
very  difficult  to  determine.  Very  little  data  exist  on  the  altitude  pro- 
files of  fogs.  Figure  18  Illustrates  typical  profiles  for  various  fogs 
as  taken  from  Jlusto  [29] . An  analytic  expression  which  approximates  the 
sea  fog  profile  is  the  following: 

8<Z)  * " zp)  U<Zc  " Z)  <62) 

where  z is  the  cut-off  for  the  fog  and  U(z  - z)  is  the  Heaviside 
c c 

function . 


[291 

J.  E.  Jlusto,  Project  Fog  Drops,  Investigation  of  Warm  Fog 
Properties  and  Fog  Modification  Concepts,  NASA  CR-72,  1964. 
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FIGURE  18.  VERTICAL  DISTRIBUTION  OF  LIQUID  WATER  CONTENT  IN  FOG 
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NATURAL  RADIATION 


In  this  section  we  shall  discuss  the  natural  sources  of  visible  and 
infrared  radiation  with  emphasis  on  nighttime  conditions. 

5.1  SOLAR  SPECTRUM 

The  extraterrestrial  value  of  the  solar  spectrum  at  the  top  of  Earth's 
atmosphere  h. s been  measured  several  times  in  the  last  few  years.  In  the 
past,  measurements  of  the  solar  spectrum  were  conducted  on  the  Earth's 
surface  for  various  air  masses  and  then  extrapolations  were  made  to  esti- 
mate the  spectrum  at  the  top  of  the  atmosphere.  A typical  set  of  solar 
values  at  sea  level  are  given  in  the  Handbook  of  Geophysics  [30] . With 
the  accessibility  of  high  flying  instrumented  aircraft  and  spacecraft, 
good  extraterrestrial  values  of  the  solar  spectrum  were  obtained  by 
Arvesen  et  al.  [31]  and  Thekaekara  and  Drummond  [32].  The  NASA  standard 
values  of  Thekaekara  and  Drummond  are  those  which  one  shall  use  in  our 
analysis.  It  is  not  known  to  what  extent  the  solar  spectrum  changes 
spectrally  over  an  extended  period  of  time.  Solar  flares  and  the  eleven 
year  solar  cycle  do  cause  changes  in  the  total  output  but  just  how  the 
spectrum  changes  is  part  of  continuing  astrophysical  research.  It  is 
generally  believed  however,  that  the  variations  in  the  visible  and  near 
infrared  parts  cf  the  solar  spectrum  are  not  greater  than  the  geometric 
variation  from  perihelion  to  aphelion.  The  latter  changes  by  about  seven 
per  cent  and  can  easily  be  calculated. 

^^USAF  Air  Research  and  Development  Command,  Handbook  of  Geo- 
physics, The  Macmillan  Co.,  New  York,  Rev.  Ed.,  1961. 

r 31 1 

1 J.  C.  Arvesen,  R.  N.  Griffin,  Jr.  and  B.  D.  Pearson,  Jr.,  Deter- 
mination of  Extraterrestrial  Solar  Spectral  Irradiance  from  a Research 
Aircraft,  Appl.  Optics,  8,  No.  11,  1969. 

f32 1 

M.  P.  Thekaekara  and  A.  J.  Drummond,  Standard  Values  for  the 
Solar  Constant  and  Its  Spectral  Components,  Nature  Phys.  Sci.,  229,  1971, 
pp . 6-9 . 
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5.2  THE  LUNAR  SPECTRUM 


5.2.1  EXTRATERRESTRIAL  VALUES 

During  the  daytime  the  sun  overwhelms  every  other  source  of  natural 
radiation.  It  has  no  intrinsic  radiative  properties  of  its  own  and  all 
radiation  which  emanates  from  the  moon  arises  from  reflected  sunlight, 
absorbed  sunlight  being  reemitted  by  a heated  lunar  surface,  or  earth- 
shine  being  reflected  by  the  moon.  Since  we  are  concerned  primarily  with 
the  visible  and  near  infrared  part  of  the  spectrum,  we  shall  only  con- 
sider the  directly  reflected  sunlight.  We  shall  now  see  how  to  calculate 
the  lunar  spectrum  at  the  top  of  Earth's  atmosphere  given  the  solar 
spectrum. 

We  can  define  the  geometrical  albedo  of  the  moon  as  the  ratio  of  the 
exitance  at  a phase  angle  of  zero  (full  moon)  to  the  total  solar  irraii- 
ance  on  the  moon,  i.e., 

p (geometrical  albedo)  ^ E(o)/Eq  (63) 

The  Bond  albedo  is  defined  as  the  ratio  of  the  radiation  reflected  over 
all  angles  to  the  solar  input,  i.e., 

2n 

A(Bond  albedo)  " ~ J"  E(a)sina  da  (64) 

o • 
o 

but,  this  can  be  written  as 


where 


(65) 


f ' 

q = 2 j $(a)sina  dot  (66) 

o 
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is  called  the  phase  integral.  The  30-called  phase  function  of  the  moon 
has  been  the  subject  of  much  investigation  over  the  years.  Early  measure- 
ments of  the  relative  intensity  of  moonlight  as  a function  of  lunar  phase 
were  made  by  Russell  [33],  Rougier  [34],  and  Bullrich  [35].  Figure  19 
illustrates  the  variation  in  a lunation  according  to  Russell.  Lane  and 
Irvine  [36]  determined  values  of  the  phase  function  on  terms  of  wavelength. 


FIGURE  19.  VARIATION  IN  A LUNATION 


t33]H.  N.  Russell,  Ap.  J.  43,  1916,  p.  117. 
f 34 1 

1 JM.  G.  Rougier,  Ann.  Obs.  Strassbourg,  2,  1933,  p.  205. 

^33^K.  Bullrich,  Ber.  Deutsch.  Wetterd.,  U.S.  Zone  No.  4,  1948. 
t^A.  P.  Lane  and  W.  M.  Irvine,  Astron.  J.  78,  No.  3,  1973,  p.  267. 
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For  our  studies  however,  this  degree  of  refinement  is  not  necessary  and 
we  shall  use  Eullrich's  values  for  all  wavelengths.  These  are  presented 
in  Table  1. 


TABLE  1.  THE  INTEGRATED  PHASE-CURVE  OF  THE  MOON 


Phase 

Angle 

Relative 

Intensity 

Phase 

Angle 

Re] ative 

Intensity 

Rougier 

Bullrich 

Rougier 

Bullrich 

0° 

100 

100 

90° 

8.0 

6.7 

10 

72.5 

73.2 

100 

5.6 

4.7 

20 

57.8 

56.0 

110 

3.9 

3.6 

30 

43.7 

42.3 

120 

2.5 

2.4 

40 

33.9 

32.0 

130 

1.6 

1.2 

50 

26.3 

23.3 

140 

0.9 

0.9 

60 

20.8 

16.7 

150 

0.4 

0.4 

70 

16.6 

12.4 

160 

0.002 

0.002 

80 

11.5 

8.7 

Values  of  the  geometric  albedo,  p of  the  moon  were  taken  from  data 
of  Condron  et  al.  [37]  and  of  Lane  and  Irvine  [36].  The  geometric  albedo 
[37]  is  shown  in  Figure  20. 

In  order  to  determine  the  extraterrestrial  value  of  the  lunar  spec- 
trum at  the  top  of  Earth's  atmosphere  we  must  know  the  mean  angle  of  the 

moon  subtended  at  Earth.  The  mean  half  angle  6 Is 

m 

= 932.7"  (67) 


[371 

T.  P.  Condron,  J.  J.  Lovett,  W.  H.  Barnes,  L.  Marcotte,  and 
R.  Nadile,  Gemini  7 Lunar  Measurements,  AFCRL-68-0438,  Sept.,  1968. 
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Thus,  the  geometric  albedo,  p is  given  by 


(68) 


where  ft  is  the  lunar  solid  angle  and  E is  the  extraterrestrial  lunar 

cm 

irradiance  at  the  top  of  Earth’s  atmosphere.  Hence,  we  have,  in  general 


Eom(A’x)  = I E0(A’x)  P(A)  P('  >v) 


(69) 


where  ft/ it  = 2.04472  x 10  The  solar  irradiance,  EQ(x»x)  depends  only 
slightly  on  phase  angle,  x and  we  shall  neglect  this  small  variation  in 
our  present  work.  The  phase  function,  P(',0  depends  little  on  the  wave- 
length X so  that  we  can  rewrite  Equation  (69)  as 

E (A,X)  = 2.04472  x 10"5  E (A)  p(A)  P(>)  (70) 

OUi  O 


where  P(o)  = 1.  Using  Equation  (70)  we  calculated  the  lunar  spectrum  at 
the  top  of  the  atmosphere  for  a full  moon.  This  is  illustrated  in 
Figure  21  along  with  the  solar  spectrum. 

5.2.2  SEA  LEVEL  VALUES  (MEASURED) 

The  sea  level  value  of  the  irradiance  due  to  the  moon  has  been 
measured  by  a number  of  investigators  over  the  years.  In  the  period 
1943-1947,  Brown  [38]  made  photometric  measurements  of  daytime  and  night- 
time illumination  as  a function  of  solar  altitude  throughout  the  world. 
Figure  22  gives  the  illuminance  as  a function  of  solar  and  lunar  altitude. 
Pleitcr  and  Morley  [39]  and  Vatsia  et  al.  [40]  measured  spectral  irradiance 

f 381 

JD.  R.  E.  Brown,  Natural  Illumination  Charts,  Report  No.  374-1, 
Dept,  of  the  Navy,  Bureau  of  Ships,  Washington,  D.C.,  1952. 
r 39  ] 

D.  Pleiter  and  G.  A.  Morley,  Night  Radiation  Measurements, 

Defense  Documentation  Center,  AD  865252. 

L.  Vatsia,  U.  K.  Stich,  and  D.  Dunlap,  J.  Opt.  Soc.  Am.  59, 

483  (1969). 
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LUMINOUS  INCIDANCE  (FOOTCANDLES) 


FIGURE  22.  NATURAL  LUMINOUS  INCIDENCE  (ILLUMINANCE)  AS  A FUNCTION 
OF  SOLAR  AND  LUNAR  ALTITUDE  AND  FOR  VARIOUS  PHASES  OF 
THE  MOON  AS  MEASURED  BY  BROWN 
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and  spectral  radiance  at  Lake  Montauban,  Canada  for  various  lunar  alti- 
tudes and  per  cent  cloud  cover.  Biberman  [41]  presents  irradiance  values 
for  moonlight  plus  airglow  for  various  lunar  phases  and  air  masses.  The 
measured  values  for  Vatsia  et  al.  [40]  are  shown  in  Figure  23.  Biberman' 8 
results  are  illustrated  in  Figures  24,  25,  26,  and  27. 

5.2.3  SEA  LEVEL  VALUES  (CALCULATED) 

Although  it  is  desirable  to  have  measurements  of  the  radiance  and 
irradiance  of  the  moon,  there  is  a practical  limit  to  the  number  of 
measurements  which  can  be  made  in  order  to  account  for  all  possible 
atmospheric  conditions . For  this  reason  it  is  necessary  to  have  an 
atmospheric-radiative-transfer  model  which  will  allow  one  to  calculate 
the  radiation  field  in  a realistic  atmosphere  for  a wide  range  of  environ- 
mental conditions.  Such  a model  was  developed  by  Turner  [42]  for  use  in 
NASA's  processing  of  multispectral  data  obtained  by  Earth  satellites. 

The  model  allows  one  to  calculate  the  spectral  irradiance  on  the  Earth’s 
surface  due  to  the  sun  or  moon  in  terms  of  solar  or  lunar  zenith  angle, 
lunar  phase,  surface  reflectance,  and  horizontal  visual  range. 

The  direct  radiation  from  the  moon  is  attenuated  exponentially  by 
the  atmosphere.  That  part  of  the  radiation  which  is  not  absorbed  is 
scattered  one  or  more  times  and  produces  what  is  called  the  diffuse  com- 
ponent. Figure  28  illustrates  the  variation  of  the  diffuse  lunar  irradi- 
ance at  the  Earth's  surface  in  terms  of  the  lunar  zenith  angle  for  a full 
moon.  The  atmospheric  states  considered  range  from  a very  hazy  atmos- 
phere of  visibility  2 km  to  the  clearest  possible  atmosphere  with  a visi- 
bility of  336  km.  Values  of  the  irradiance  for  angles  greater  than  ^75° 
are  not  tc  be  considered  as  being  extremely  reliable  since  refraction  of 

— — 

L.  M.  Biberman,  Natural  Levels  of  Illumination  and  Irradiance, 
Photoelectronic  Imaging  Devices,  Vol.  1,  Chapt.  3,  Plenum  Press,  N.Y., 
1971. 

T421 

1 JR.  E.  Turner,  Radiative-Transfer  in  Real  Atmospheres,  ERIM 
190100-24-T,  Environmental  Research  Institute  of  Michigan,  1974. 
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FIGURE  23.  NIGHT  SKY  RADIANT  STERANCE  SPECTRA  FOR  VARIOUS  ALTITUDES 
OF  THE  MOON  AS  MEASURED  BY  VATSIA,  STICH  AND  DUNLAP 
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FIGURE  24.  IRRADIANCE  FROM  FULL  MOONLIGHT  PLUS  AIRGLOW 
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FIGURE  25.  IRRADIANCE  FROM  0,3  x FULL  MOONLI'JHT  PLUS  AIRGLOW 


60 


J2 


Wavelength,  nanometers 


FIGURE  26.  IRRADIANCE  FROM  0.1  x FULL  MOONLIGHT  PLUS  AIRGLOW 
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FIGURE  27.  IRRADIANCE  FROM  0.03  x FULL  MOONLIGHT  PLUS  AIRGLOW 
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Diffuse  Irradlance  x 10J  (milliwatts  per  cm  per  micrometer) 


rr-T-i-j 


Jpj 

r 

" r 

• 

: 

j.  - - 

1 

-rr: 

r 

7-| ! r 

7J 

s 

....: .. 

* 

.... 

i 

. ..j.. .... 

: 

| _ 

i n 

..  .... 

FIGURE 

28 

...r. 

• ; i 

; ' 

!‘  - 

).  V-2 

■-*.  • 

...  \ ! 

km  ...  L — 

j .. 

...._. 

■ 

i.... 



: 

1 

...  . 

DIFFUSE  LUNAR  IRRADIANCE  AT  THE  EARTH’S 
SURFACE  VS.  LUNAR  ZENITH  ANGLE 

V “ Visibility 

i 

i 

1 

•• 

i ; 

....: 

? 

Wavelength  = 

.55 

licrometers 

J 

: 

. 

; 

. ..; ... 

; 

; 

v. 

: 

..  — 

N< 

'll  ' 

Hazy  Atmospnere 
Background  Reflectance 

= 0.0 

'"7~ 
^ . 

. \ : 
: 

L 

r • 

i' 

Lunar  Phase  Angle 

= 0.0 

. 

L .. 

• 

■ 

' \ 1 

r ■ 

j 

1 

; ■ 

: j 

• -i 

rr 

! 

i 

1 

> 

1 

' 

- . 

I i j 

■ \ . 

L;  .. 

j 

■ i 

i • 

• 

\ 

r , 

1 

1 

i 

’ i j 

V ' 

r 

i 

j 

\ 

7 ' 

| 

*j 

. ; 

i : 

■ : 

| 

j 

1 

• - 

; 

\ 

r ] 

i 

: 

;... 

! . i 

! 

! 

: 

■-  i- 

\ 

! 

: 

. ; j 1 

j 

f: 

V-  • 

: 

1 

" !*" 

— : — 

l_h 

- -- 

. -i  • 

! - 

. 

: 

. • 

\ 

. 1 .. 

1 — 

. 

rn 

_ • ^ 

. . .. 

i : 

1 1 1 

Li- 

L 1.. 

“ * 

] 

v? 

! - 

r* 

• 

. ..! 

r"; 

j 

\ : 

: i i . 

!• 

. 1 

! 

; 

i 

:: 

. - i 

! S 

s : 

. j _ 

; i • | 

! 

i 

* : 

1 

i 

”71 

• 7' 

rJ.J 

; i I 

fe 

t 

• 

..J._ 

. . 

! 

_ .1 

■ 

. 

“T“ 

; 

"T  • 
..._ 

* ■ 

■‘7 
. • . 
i 
* 

i 

" j 
1 

t 

■ ! 1 

_ 

□ . 
■ 

7 

•— j- 

! 

■ 

err;... 
: ! 

** 

: 

r ^ 
_ 

X 

i 

\ 

: ; . 

; 

u -i  ! ’ 

V-23  k 

. 

1 

. 

J, 

i_ 

' 

■ I’  • 
; 

I 

......  . 

: 

..  . 

i 

' 

1 

n 

i 

’ 

n 

j » 

i 1 

LU— L:  .. 

.!.. 

:: 

..  — - 

i 

ITL 

H 

. _ L. 
! 

.... 

_A 

■ , 
....... 

......  j 1 

' ’ 

j 

T 

. 

: 

: : 

... 

j . 

.....  . 

: i : 

....:  i . 

• i ! 

..... 

.... 

1 

! ' 

• 

. : . 

... 

rr 

! 

‘ 

rj 



\ 

r ; 

V-326  km 

! ; 

. ! . 
_.!  :... 

. ; 

; 

L, .. 

• 

— 

. 

: 

• -■ 
- 

. 

V \ 

V ' 

0 

... 

. 

7]'i“ 

; 

i . 

..  i 

• t 

! * r — 

i 

. _.L- 

■ 

* 

... 

. 

the  atmosphere  was  not  used  in  the  modal.  Figure  29  illustrates  the  same 
variation  except  cor  a quarter  moon. 

The  total  lunar  irradiance  is  the  sum  of  the  diffuse  component  and 
the  directly  attenuated  radiation.  Figures  3Q  and  31  illustrate  the 
variation  of  total  lunar  irradiance  with  lunar  zenith  angle  for  the  full 
moon  and  the  quarter  moon.  As  can  be  seen,  the  dependence  on  visibility 
is  small.  This  is  a result  of  the  sum  of  the  direct  and  diffuse  compon- 
ents. The  diffuse  radiation  is  greater  for  hazy  atmospheres  and  the 
direct  component  is  less  than  for  clear  atmospheres.  Ilence,  the  sum  of 
the  two  opposing  effects  (attenuation  and  scattering)  combine  in  such  a 
way  as  to  minimize  the  variation  with  visibility.  This  effect  is  depicted 
in  Figures  32  and  33  for  the  full  and  quarter  moons. 

Transmittance  is  sometimes  a more  meaningful  quantity  than  visi- 
bility. Let  us  consider  the  variation  of  the  diffuse  and  total  irradi- 
ance with  transmittance  at  a wavelength  corresponding  to  the  peak  value 
of  the  visual  response  of  the  human  eye.  Figure  34  illustrates  this 
variaticn.  For  extremely  thick  atmospheres,  as  in  fogs,  the  total  and 
diffuse  components  have  the  same  value  for  a given  surface.  When  the 
transmittance  reaches  a value  of  about  one  per  cent  however,  the  diffuse 
irradiance  decreases  sharply  with  a further  increase  in  transmittance. 

At  a transmittance  of  one  all  the  total  irradiance  curves  have  the  same 
value,  which  are  independent  of  reflectance  as  they  should  be  since  the 
total  irradiance  becomes  equal  to  the  direct  irradiance  at  the  top  of  the 
atmosphere.  It  is  interesting  that  the  total  irradiance  and  the  diffuse 
irradiance  can  be  greater  than  the  input  value  at  the  top  of  the  atmos- 
phere. This  arises  because  for  highly  reflecting  surfaces  the  diffuse 
radiation  reflects  back  and  forth  between  the  surface  and  the  overlying 
turbid  atmosphere  and  the  radiation  is  effectively  trapped. 

The  dependence  on  reflectance  can  be  understood  more  easily  by 
looking  at  Figure  35.  For  atmospheres  which  are  not  extremely  thick 
there  is  almost  a linear  relationship  connecting  the  diffuse  irradiance 
with  reflectance  but  for  optically  thick  atmospheres  the  relationship  is 
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FICURE  29 


DIFFUSE  LUNAR  IRRADIANCE  AT  THE  EARTH'S 
SURFACE  VS.  LUNAR  ZENITH  ANGLE 
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non-linear.  It  is  interesting  to  note  that  at  least  in  this  particular 
case  the  diffuse  irradiance  is  independent  of  transmittance  for  a 
reflectance  near  0.70. 

Finally,  we  can  investigate  the  dependence  of  total  irradiance  on 
reflectance  by  analyzing  Figure  36.  Here  we  see  that  for  ordinary  atmos- 
pheric states  with  reasonable  transmlttances  the  total  irradiance  is 
essentially  independent  of  reflectance.  For  foggy  atmospheres,  however, 
there  is  a distinct  non-linear  relationship.  For  low  values  of  reflec- 
tance and  as  the  atmosphere  becomes  clear  more  radiation  is  able  to  reach 
the  surface  and  therefore  the  total  irradiance  increases.  On  the  other 
hand,  for  high  reflectance  values,  and  as  the  atmosphere  becomes  clear, 
more  radiation  is  able  to  leave  the  atmosphere  by  reflection  and  hence 
the  total  irradiance  decreases . As  in  the  case  of  diffuse  Irradiance 
there  exists  a small  region  of  reflectance  values  such  that  the  irradi- 
ance is  essentially  independent  of  transmittance. 

In  this  section  we  have  seen  how  the  lunar  irradiance  varies  with 
the  normal  environmental  parameters.  The  spectral  dependence  of  irradi- 
ance is  that  of  the  extraterrestrial  lunar  value  modified  by  the  basic 
spectral  dependence  of  whatever  atmosphere  and  surface  reflectance  exists. 
To  consider  all  of  these  functional  relationships  is  beyond  the  immediate 
scope  of  this  study.  These  investigations  are  included  in  a report  by 
Malila  et  al.  [43] . 


^ W.  A.  Malila,  R.  B.  Crane,  and  R.  E.  Turner,  Information 
Extraction  Techniques  for  Multispectral  Scanner  Data,  Report  No.  31650- 
74-T,  Willow  Run  Laboratories,  Ann  Arbor,  1972. 
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FIGURE  36 


TOTAL  LUNAR  IRRADIANCE  AT  THE  EARTH’S 
SURFACE  VS.  BACKGROUND  REFLECTANCE 

T “ Transmittance 
o 

Wavelength  ■ .55  Micrometers 
Hazy  Atmosphere 
Lunar  Zenith  Angle  “ 30.0 
Lunar  Phase  Angle  *>  0.0 


6 

AFTIFICIAL  SOURCES 


In  this  section  we  shall  consider  the  radiation  fielcl  in  an  atmos- 
phere composed  of  isotropic  point  sources  (flares) . The  general  mathe- 
matical problem  of  determining  the  complete  multiply  scattered  field  is 
quite  complicated.  Instead,  we  shall  concentrate  on  the  most  Important 
terms,  the  directly  attenuated  irradiance  on  the  horizontal  surface  and 
singly-scattered  irradiance. 

6.1  POINT  SOURCE 

The  general  three-dimensional  radiative-transfer  equation  for  a 
point  source  is 

p(O.Q')  L(r,6')dQ'  = 6(r-r  ) (71) 

4tr  o 

fi 

where  O'  is  the  direction  of  the  incoming  ray  and  0 is  the  direction  of 
the  received  ray.  S is  the  spectral  power  output  of  the  source  in  watts 
and  L(r,0)  is  the  spectral  radiance  at  point  r.  The  spectral  irradiance 
on  the  surface  is  given  by 


Q . VL(r ,0)  + tc(r)  L(r,fl) 


-Mil  f 
An  J 


E(R) 


j n*«'  Klt.O^dO' 
O' 


(72) 


where  R is  a vector  fro r.  the  origin  of  coordinates  to  the  point  in 
question  located  on  the  surface.  The  vector  n is  the  surface  normal 
vector . 

Actually,  Equation  (71)  is  used  to  find  the  Green's  function  for  the 
radiative-transfer  equation.  This  function  is  then  the  solution  of  the 
basic  equation  for  a point  source. 


73 


6.2  DISTRIBUTED  SOURCE 

The  solution  for  a spatially  distributed  source  is  given  by 


L(r,0) 


1 1 p(v 


V 


G(f,G,1? ,n  )dr  dfi 
0 0 o o 


(73) 


where  P(?o.0o)  i8  the  spatial  distribution  of  the  source  and  G(f,n,t  ,(2  ) 
is  the  Green's  function.  Thus,  we  cun  integrate  over  the  source  coordi- 
nates and  obtain  the  radiance,  L(F,0)  at  any  point  ?.  Examples  of  dis- 
tributed sources  are:  clouds,  buildings,  mountains,  burning  vehicles, 

and  other  fires. 

6.3  THE  RADIATION  FIELD 

We  shall  now  consider  the  radiation  field  from  a point  source.  The 

geometry  is  illustrated  in  Figure  37.  The  radiative-transfer  equation 
is 

fi.VL  + <(z)  L(x,y,z,u ,4>) 


2tt  1 

f I P(u,<t>,v',<p’)  L(x,y,z,y'5(f')du'd^'  - 6(x)6(y)6(z-z  ) 

o -1 

(74) 


We  separate  the  radiation  fields  into  two  parts,  the  directly  attenuated 
and  the  singly-scattered  terms,  i.e., 

L(x*y*z»P  *$)  ■ L<j(x,y,z,y  ,<p)  + Lg(x,y,z,u,<j>)  ( 


Taking  the  Fourier  transform  of  Equation  (74)  and  performing  various 
mathematical  operations  we  get  fcr  the  direct  radiance 


T-*d(x,y, *,**,♦) 


S -(t(zq)-t(z)]/p 

4iry  e 


6[X  - y <2-Z0)]  6[y  (Z_Z0)]  U<2-20)  (76) 
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for  the  upward  directed  radiance  and 


s -[t(z)-t(zo)]/p 

L^.y.z.-H.l)  '^e  x 

Six  - (zQ-z)]  5[y  “ “ (zq-z)]  U(z0-z)  (77) 


for  the  downward  radiance  where  t(z)  is  the  optical  depth  at  altitude  z 
and  £ and  u are  direction  cosines  with  respect  to  the  x and  y axes. 
Integrating  the  downward  radiance  over  a hemisphere  gives  us  the  down- 
ward attenuated  irradiance  due  to  a point  source. 


VR> 


z S e 
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3/2 


(78) 


where  R is  the  distance  from  the  origin  along  the  surface. 

Now,  integrating  the  radiative-transfer  equation  once  and  performing 
considerable  mathematical  analysis  leads  to  the  following  solution  for 
the  radiance 


L (R,n,0,-y) 

o 


-[Tq-T (** ) ] /p  - 


16n2p  J 


p (cos  X)e 


t(zq)-t(z') 


z'-z. 


O -J 


B(z')dz' 
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where 


\/r2  + z2  - 2z  z'  + ~ + 
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2 ■ 2Rz'  — cos  n 


(80) 
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and 


cos  x 


cos  n 


(81) 


The  spectral  irradiance  due  to  single  scattering  from  a point  source 
is  then  given  by 


2*  1 

Es<*>  “|  I u Lg(R,n,  0,-y)dydn  (82) 

o o 


6.4  RESULTS  OF  CALCULATIONS 

We  are  now  going  to  see  the  effect  of  one  or  more  flares  on  the 
illumination  of  the  ground.  Also,  we  shall  illustrate  the  illumination 
due  to  the  moon. 

A typical  scene  is  depicted  in  Figure  38.  We  have  moonlight  and 


FIGURE  38.  NATURAL  AND  ARTIFICIAL  ILLUMINATION  CONFIGURATION 
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artificial  flares  illuminating  the  surface.  The  moonlight  is  the  result 
of  an  extraterrestrial  lunar  irradiance  being  scattered  by  the  atmosphere. 
The  scattering  is  multiple  scattering  and  we  also  include  the  effect  of 
reflection  from  the  surface.  For  the  flares  we  calculate  the  directly 
attenuated  and  the  singly  scattered  components. 

First,  we  should  consider  a single  flare  in  a hazy  atmosphere  at 
three  different  altitudes.  The  total  (direct  plus  singly-scattered) 
irradiance  is  illustrated  in  Figure  39.  The  results  are  normalized  to  a 
spectral  power  output  of  one  watt  per  micrometer.  Any  other  flare  output 
is  just  multiplied  by  these  values.  It  is  Interesting  to  note  that  for 
the  higher  flares  the  irradiance  is  less  than  for  lower  flares  when  one 
is  near  the  origin  (directly  beneath  the  flare).  At  greater  distances 
from  the  origin,  however,  it  is  reversed.  This  is  due  to  tne  fact  that 
the  lower  flares  are  more  deeply  embedded  in  :he  heavy  haze  layer  near 
the  surface  where  more  scattering  occurs. 

He  now  look  at  the  case  of  a uniform  fog  which  has  a constant  pro- 
file up  to  an  altitude  of  0.1  km.  The  visibility  is  only  100  meters. 

This  case  is  Illustrated  in  Figure  40.  Here  we  see  the  same  effect  as  in 
the  haze  except  that  it  is  more  severe,  i.e.,  the  lower  sources  provide 
more  illumination  at  greater  distances  than  the  higher  sources.  An  excep- 
tion is  to  be  noted  however,  for  the  very  low  source  of  0.05  km.  The 
other  three  sources  are  above  the  fog  and  scatter  the  radiation  more, 
whereas  the  lowest  source  lies  within  the  fog  and  the  radiation  is  attenu- 
ated strongly. 

Figure  41  depicts  the  case  of  a much  heavier  fog.  Here  we  see  the 
same  effect  as  before  only  more  so.  The  visibility  in  this  case  is 
10  meters. 

He  can  now  look  at  the  effect  of  an  altitude  variation  of  the  source. 
Considering  the  100  meter  visibility  fog  the  effect  of  this  variation  is 
shown  in  Figure  42.  It  should  be  noted  that  when  we  are  at  the  origin 
the  irradiance  decreases  monotonically  as  the  source  moves  farther  away. 
If,  however,  we  are  at  greater  distances  along  the  surface  then  Lhe 
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SINGLY  SCATTERED  I RRAD LANCE  (milliwatts  per  cm  per  micrometer) 


FIGURE  39 


DIRECT  PLUS  SINGLY  SCATTERED  IRRADIANCE 
GENERATED  BY  A POINT  SOURCE  OF  RADIATION 
VS.  DISTANCE  ON  THE  GROUND  FROM  THE  SOURCE 

Source  Strenpch  * 1 watt  per  micrometer 

Z#  “ Altitude  of  source 

Wavelength  ■ .55  micrometers 

Hazy  Atmosphere 

Visibility  ■ 2 km 
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irradiance  decreases  down  to  a minimum  value,  then  increases  to  a maximum 
and  finally  decreases  monotonically  thereafter.  Here  we  are  seeing  the 
trade-off  between  attenuation  and  scattering.  For  low  altitude  sources 
the  radiation  is  primarily  being  attenuated,  hence  the  decrease.  As  the 
source  height  increases  the  singly-scattered  radiation  is  more  important 
and  the  irradiance  begins  to  increase.  There  is  a limit  to  its  increase, 
of  course,  since  eventually  we  have  a source  well  above  the  fog  layer  and 
the  amount  of  radiation  decreases  as  the  inverse  square  of  the  distance 
from  the  source.  Therefore,  a peak  value  is  reached  and  a monotonic 
decrease  proceeds  from  that  point  on. 

We  shall  now  consider  the  case  of  a fixed  point  source  at  an  alti- 
tude of  200  meters  embedded  in  a fog  with  varying  visibility.  This  is 
illustrated  in  Figure  43.  For  very  low  visibilities  the  attenuation  is 
quite  strong  and  almost  no  radiation  arrives  at  the  surface.  As  the 
visibility  increases  there  is  less  attenuation  and  hence  a sharp  increase 
in  irradiance  occurs.  However,  as  the  visibility  continue^  to  increase 
the  scattered  radiation  decreases  and  therefore  the  total  irradiance 
decreases . 

Figure  44  illustrates  the  same  effect  for  a distance  of  1 km  from 
the  source.  The  same  basic  variation  occurs  as  in  the  previous  case 
except  the  total  magnitude  is  less . 

We  should  investigate  the  spectral  dependence  of  the  irradiance. 

We  consider  a spectrally  flat  (white)  flare  first  in  order  to  see  the 
spectral  dependence  of  the  atmosphere.  Figure  45  illustrates  this  varia- 
tion for  a hazy  atmosphere  with  a visibility  of  13  km.  There  is  little 
spectral  variation  in  the  directly  attenuated  irradiance  but  the  singly 
scattered  component  does  decrease  monotonically  with  increasing  wave- 
length . 

Figure  46  illustrates  the  same  dependence  except  for  a distance  of 
one  kilometer  from  the  origin.  Since  we  are  farther  from  the  source  the 
attenuation  and  scattering  is  greater,  hence  the  greater  degree  of  spec- 
tral dependence. 
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SINGLY  SCATTERED  IRRADX^JCE  (milliwatts  par  cm  par  mlcromatar) 


FIGURE  44 


i . I l ...  I 

DIRECT  PLUS  SINGLY  SCATTERED  IRRADIANCE 
GENERATED  BY  A POINT  SOURCE  OF  RADIATION 
VS.  VISIBILITY  IN  A FOG 

Source  Strength  ■ 1 watt  per  micrometer 
Altitude  of  source  - .2  km 
Distance  on  '-he  ground  from  the  source-1 . Okm 
Wavelength  “ .55  micrometers 

Homogeneous  Fog  Confined  Below  .1  km 
Altitude 


Figures  47  and  48  demonstrate  similar  spectral  variations  for  fogs 
instead  of  hazes.  Since  fog  droplets  are  much  larger  than  the  particu- 
lates in  a hazt  the  spectral  dependence  is  much  less  as  we  saw  in 
Section  4. 

We  can  now  consider  the  spectral  variation  for  actual  flares.  Data 
were  provided  us  on  actual  flare  spectra  by  Picatinny  Arsenal.  The 
spectra  are  presented  in  Figures  49  and  50.  The  change  in  the  actual 
flare  spectrum  for  Mark  45  flare  is  illustrated  in  Figure  51  for  a hazy 
atmosphere.  The  direct  and  singly-scattered  irradiance  at  first  appear 
to  be  the  same.  Closer  inspection  indicates,  however,  that  the  direct 
irradiance  is  slightly  greater  in  the  infrared  relative  to  the  singly- 
scattered  irradiance. 

Summarizing  these  results,  we  can  say  that  for  "clean"  fogs  in  which 
there  are  no  absorbing  contaminants  with  an  intrinsic  spectral  dependence, 
then  the  spectral  variation  in  Irradiance  is  very  small.  For  hazes  or 
fogs  contaminated  with  dust  particles  the  spectral  variation  could  be 
quite  different.  These  atmospheric  states  should  be  Investigated. 

Next,  we  shall  investigate  multiple  flares  in  a realistic  nighttime 
environment . 

First,  we  shall  consider  three  flares  in  a right  triangle  configu- 
ration as  illustrated  in  Figure  52.  They  are  at  altitudes  of  0.2  km, 

0.1  km,  and  0.4  km  in  a very  hazy  atmosphere  with  a visibility  of  2 km. 

It  should  be  noted  that  the  highest  flare  has  a high  contour  closer  to 
the  center  than  for  the  other  flares.  This  effect  was  explained  in  the 
description  of  the  single-flare  graphs. 

Figure  53  illustrates  the  contours  for  the  same  situation  as  in 
Figure  52  except  that  the  flares  are  twice  as  far  apart. 

We  can  now  look  at  the  case  of  a fog.  Figure  54  shows  three  flares 
in  a foggy  atmosphere.  The  fog  has  a constant  profile  up  to  a height  of 
100  meters.  One  flare  is  at  an  altitude  of  50  m,  another  at  the  top  of 
the  fog  layer,  and  the  third  lies  50  meters  above  the  fog  layer.  Here  we 
see  a very  interesting  effect.  The  high-altitude  flare  produces  a wide 
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IRJRADIANCE  (milliwatt*  par  cm  par  mlcromatar) 


~ I j~^~:  w'WUS  48  “ • - | — -rj-  - ~ j - 

DIRECT  AND  SINGLY  SCATTERED  IRRADIANCE 
GENERATED  BY  A POINT  SOURCE  OF  RADIATION 
VS.  WAVELENGTH 

Source  Strength  - 1 watt  per  micrometer 
(for  all  wavelengths) 

Altitude  of  the  source  ■ .2  ka 
Distance  on  the  ground  fro*  the  source*  l.Oka 
Homogeneous  Fog  Confined  Below  .1  km  Altitude 
Visibility  (in  fog)  » .1  ka 


Singly  Scattered 
Irradiance 


NOTE  CHANCE  IN  SCALE 


Direct  Irradiance 


WAVELENGTH  (MICROMETERS) 


per  micrometer) 


FIGURE  58 


POWER  OUTPUT  OF  SR  485  FLARE 
EXTRAPOLATED  TO  MARK  45  SIZE 
AMD  BURNING  RATE  VS.  WAVELENGTH 
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FIGURE  SI 


DIRECT  AND  SINGLY  SCATTERED  IRRADiAi.CE 
GENERATED  BY  A MARK  AS  FLARE  VS.  WAVELENGTH 

Altitude  of  Flare  - .2  kn 

Distance  on  the  ground  from  the  flare-1.0  km 

Hazy  Atmosphere 

Visibility  - 13.0  km 


DIRECT  IRRAD  LANCE 


WAVELENGTH  (MICROMETERS) 


Wavelength 
Hazy  Atraosi 
Visibility 

Altitude  of  Fla: 

.2 

.1  km  .U 


1*5  FLARES 


eters 

d Belov  .1  km 


L km 


Altitude  of  Flares  = .10  km 
.05  km  .15  km 


distribution  of  radiation  over  an  extended  area  whereas  the  lower  altitude 
flares  have  a much  more  concentrated  distribution.  This  is  due  to  the 
fact  that  scattering  is  a more  important  effect  than  attenuation  for  the 
higher  altitude  case. 

Figure  55  illustrates  the  distribution  due  to  four  flares  in  a fog. 

As  in  the  previous  case,  the  two  flares  above  the  fog  layer  produce  a more 
extended  distribution  than  the  two  flares  below  the  fog  layer. 

Finally,  in  Figure  56  we  present  four  flares  in  ar.  extremely  dense 
fog  with  a visibility  of  only  10  meters.  The  height  of  the  fog  is  30  m 
and  two  of  the  flares  are  within  the  fog  whereas  two  are  above.  Again, 
we  see  the  same  effect  as  before.  The  lowest  flare  decreases  most  rapidly. 

The  flares  are  normalized  to  one  watt  per  micrometer  and  the  values 
of  natural  irradlance  are  given  in  Section  5.  Hence  any  flare  output  can 
be  compared  with  any  natural  background. 

6.5  ILLUMINATION  MODEL 

A complete  model  for  the  determination  of  the  illumination  of  a 
scene  as  a result  of  natural  and  artificial  radiation  sources  can  be 
developed  based  upon  the  analysis  in  Sections  5 and  6.  Figure  57  is 
meant  to  illustrate  a more  complete  model  for  surface  illumination.  It 
should  not  be  Interpreted  as  a single  computer  program.  The  present 
capability  is  as  follows:  The  scene  can  oe  defined  in  terms  of  surface 
reflectance  for  a uniform,  flat  terrain.  Any  number  o;  point  sources 
may  be  used.  Meteorological  data  involving  atmospheric  pressure,  hori- 
zontal visibility,  and  fog  height  can  be  specified  to  determine  the 
atmospheric  state.  The  date  and  time  should  be  specified  in  order  to 
determine  the  extraterrestrial  solar  or  lunar  irradlance  at  the  top  of 
the  atmosphere  and  the  solar  or  lunar  zenith  angle.  Knowing  the  wave- 
length and  the  appropriate  meteorological  data,  the  optical  parameters 
such  as  attenuation  coefficients  and  optical  depth  can  be  calculated  and 
these  parameters  are  then  used  in  the  point  source  or  natuval  source  model 
to  calculate  the  natural  and  artificial  radiation  fields. 
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FIGURE  56 


CONTOURS  OF  CONSTANT  DIRECT  PLUS  SINGLY  SCATTERED  IRRADIANCE 
GENERATED  BY  FOUR  MARK  1*5  FLARES 

Wavelength  * .55  micrometers  Altitude  of  Flares  * .015  tan  .045  kn 

Homogeneous  Fog  Confined  Belov  .03  km 

Altitude  *030  kra  .CoO  km 

Visibility  (in  fog)  * .01  km 


-3  2 

Contour  levels  are  in  units  of  10  milliwatts  per  cm  per  micrometer 
Scale  * .05  km  per  inch 
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ILLUMINATION  MODEL 


FIGURE  57.  ILLUMINATION  MODEL  FOR  NATURAL  AND  ARTIFICIAL  SOURCES 
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The  present  node Is  can  be  generalized  to  include  gaseous  absorption 
and  multiple  scattering  as  veil  as  the  effects  of  reflection  from  non- 
Lambertian  surfaces.  In  addition,  further  modeling  of  cloud  layers 
should  be  conducted  in  order  to  understand  how  the  illumination  on  the 
surface  changes  as  flares  penetrate  the  cloud  layer.  Also,  the  basic 
mathematical  procedures  have  already  been  developed  which  will  allow  one 
to  calculate  the  illumination  due  to  distributed  (non-point)  sources 
such  as  fires,  city  lights,  mountains,  and  other  natural  or  man-made 
sources . 


7 

COMPUTATIONAL  METHODS 


The  mathematical  model  used  to  calculate  the  direct  and  singly 
scattered  irradiance  generated  by  an  isotropic  point  source  is  in 
principle  an  exact  formulation.  However,  some  approximations  and  simpli- 
fications are  required  to  adapt  the  solution  of  the  3-fold  Integral  equa- 
tion for  the  singly  scattered  irradiance  to  the  computer.  This  equation 
is  written  as  follows: 


Sbt 
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(33) 


S is  the  source  strength  in  watts  per  micrometer,  is  the  single 

scattering  albedo,  k(z)  is  the  vertical  attenuation  coefficient,  p(cos  x) 

is  the  scattering  phase  function,  t(z)  Is  the  optical  depth,  u is  the 

cosine  of  the  angle  between  a vertical  line  and  a line  drawn  between  the 

scattering  center  and  the  irradiated  spot  on  the  ground  (at  a distance  R 

from  a point  on  the  ground  directly  below  the  source),  z is  the  altitude 

s 

of  the  source,  D is  the  distance  in  kilometers  from  the  source  to  the 
scattering  center,  and  u'  are  angular  variables  of  integration  for  the 
spherical  coordinate  system  described  in  Section  7.2,  and  the  factor  of 

107  in  t;  a denominator  converts  the  units  of  the  result  to  milliwatts 
per  cn^  per  micrometer. 
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7.1  SIMPLIFICATIONS  AND  APPROXIMATIONS  WHICH  REDUCE  COMPUTATIONAL  EXPENSE 

The  simplifications  and  approximations  which  are  used  in  the  com- 
puter solution  for  E fall  into  two  major  categories.  The  first  of 

ss 

these  categories  is  concerned  with  making  the  evaluation  of  the  equation 
economical.  In  particular,  the  aerosol  or  fog  scattering  phase  functions 
used  in  the  computations  are  determined  by  linear  interpolation  from 
values  in  a table.  This  result  is  then  combined  with  the  Rayleigh  phase 
function  (which  is  calculated  analytically) , with  each  phase  function 
being  weighted  by  the  ratio  of  its  scattering  coefficient  to  the  total 
scattering  coefficient.  This  interpolation  and  combination  technique 
avoids  the  very  costly  alternative  of  calculating  the  phase  function  from 
a polynomial  series  (which,  since  it  requires  a great  many  terms  to  be 
accurate,  has  been  done  just  once  to  generate  the  interpolation  table). 

The  interpolation  is  done  in  two  variables:  scattering  phase  angle  and 

wavelength.  The  Interpolation  with  respect  to  wavelength  is  done  only 
once  at  the  beginning  of  the  computations,  and  this  result  is  then  inter- 
polated with  respect  to  the  scattering  phase  angle  as  required. 

The  >"vrtical  attenuation  coefficient,  <(z),  and  the  optical  depth, 
t(z),  are  actually  functions  of  wavelength  and  visibility  as  well  as 
altitude.  For  fogs  k(z)  and  t(z)  are  computed  analytically  a;?  functions 
of  visibility  and  altitude,  but  their  dependence  on  wavelength  is  deter- 
mined by  Interpolation  from  a table,  representing  either  a selective  fog 
or  a stable  fog,  as  specified  by  an  input  to  the  computer  program.  The 
analytical  dependence  of  k(z)  and  x(z)  on  altitude  is  selected  among  four 
alternatives  according  to  another  program  input.  For  hazes  (aerosol 
scattering)  k(z)  and  x(z)  are  computed  from  a piece-wise  analytical  repre- 
sentation of  Elterman's  model.  The  dependence  of  k(z)  and  t(z)  on  wave- 
length is  first  calculated  using  an  8^  order  polynomial.  For  altitude.’, 
less  than  5 km  or  greater  than  26  km  an  exponential  variation  with  alti- 
tude is  assumed,  while  for  altitudes  between  5 km  and  26  km  another  8^ 
order  polynomial  is  used  to  represent  the  altitude  dependence  of  the  model. 
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This  altitude  dependent  part  is  then  multiplied  by  the  wavelength 
dependent  part  to  obtain  k(z)  and  t(z).  The  visibility  enters  the 
calculations  (in  a straightforward  analytical  way)  only  for  altitudes 
less  than  5 km.  The  values  of  <(z)  and  t(z)  determined  above  are  added 
to  the  values  of  k(z)  and  t(z)  for  Rayleigh  scattering,  which  are  cal- 
culated analytically.  For  hazes  the  effects  of  ozone  absorption  are 
also  included  as  determined  by  linear  interpolation  from  tables  in  wave- 
length and  altitude. 

It  is  desirable  to  perform  the  integration  with  respect  to  z 
(altitude)  as  the  innermost  integral  (to  determine  radiance  as  a function 
of  <j>'  and  u'  as  well  as  to  minimize  the  number  of  calculations  in  the 
innermost  integral),  yet  it  is  not  desirable  to  recompute  k(z)  and  t(z) 
during  each  cycle  through  this  integral.  Hence,  as  much  as  possible  the 
expressions  in  the  equation  which  involve  k(z),  t(z),  z,  and  constant 
factors  are  calculated  once  at  the  beginning  of  the  computer  program  and 
are  saved  in  an  array  from  which  they  may  be  readily  obtained  when  they 
are  needed. 

The  variables  p and  w'  are  related  as  follows: 


v 


P ' cosa  + 


sina  cos<J> ' 


(84) 


while 


D = 


(85) 


and 


cos  x 


(86) 


where  P is  the  distance  between  the  source  and  the  irradiated  spot  on  the 
ground,  a is  the  angle  between  a vertical  line  and  a line  drawn  between 
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the  source  and  this  irradiated  spot,  z is  the  altitude  of  the  scattering 
center,  and  x is  the  scattering  phase  angle.  (Note:  The  limits  of  the 

integration  in  4> * and  y'  must  correspond  to  integration  over  the  hemi- 
sphere specified  by  0 <.  p <L  1.  This  is  accomplished  by  changing  the 
signs  of  both  p and  p*  whenever  the  calculation  of  p from  <j>'  and  p' 
yields  a result  less  than  zero.  D and  cosx  are  then  calculated  using 
these  new  values  of  p and  p'.  This  corresponds  to  a reflection  of  the 
coordinates  <J>'  and  p'  through  the  origin.)  During  the  integral  with 
respect  to  <J> ' , cos<J> 1 has  to  be  recomputed  for  each  change  in  <p ' . Letting 
the  increment  in  <j>’  be  represented  by  A<J>  * , we  use  the  trigonometric 
identities  for  the  sine  and  cosine  of  the  sum  of  two  angles  to  compute 
successive  values  for  sin4>'  and  cos^'.  This  requires  that  we  compute 

sinA$*  and  cosA^'  once  and  then  take  appropriate  sums  and  differences  of 

products  between  sin<t>*  or  cos<|) ' and  sinA^1  or  cosA^': 

cos(<j>’+A<(>  ’)  = cos<|>'  cosA4>'  - sin4>'  sinA((!,  (87) 

sin <<])  *+A4) * ) = sin<t>'  cosA<J>'  + costj)*  sinA<J>'  (88) 

Although  it  appears  that  one  only  needs  to  increment  cos<J>',  one  must  also 
increment  sin<J>'  in  order  to  keep  this  sequence  going.  This  method  for 
incrementing  cos^'  is  more  efficient  than  calculating  cos<J>'  from  each 
new  value  of  <J>  * . 


Since  the  calculation  of  E has  Dilateral  symmetry,  one  need  only 

s s 


integrate  for  0 <.  4> * <.  * and  multiply  this  result  by  2. 
further  simplification  is  also  possible.  Since 

However,  a 

cos(j  - <p'  ) = sin<t>' 

(89) 

cos (^  + <j>')  * -sin^' 

(90) 

COS (TT  - <j>')  * -COSl}l, 

(91) 

10b 


the  <t>'  integration  can  be  split  into  four  more  parts . Hence,  we  inte- 
grate over  the  range  0 <_  <_  , but  add  together  four  separate  results 

for  each  value  of  <j> ' — the  straightforward  result  corresponding  to  the 
equation  for  u as  a function  of  and  u*  as  written,  and  three  other 
results  corresponding  to  replacing  cos<{> ' ii  the  equation  for  p by  sin*’, 
-sin<t>',  and  -cos^',  respectively.  The  appropriate  substitute  values  for 
cos<t>'  are  all  readily  available  from  the  incrementing  scheme  described 
above . 

7.2  SIMPLIFICATIONS  AND  APPROXIMATIONS  WHICH  MAKE  THE  EQUATION  SOLUBLE 

The  second  major  category  of  simplifications  and  approximations  used 

in  the  calculation  of  E on  the  computer  is  concerned  with  making  evalu- 

s s 

ation  of  the  equation  possible.  Since  the  integrals  cannot  be  evaluated 
analytically  in  a closed  form,  some  sort  of  numerical  integration  method 
is  required.  Two  basic  considerations  are  of  importance  here:  (1)  the 

choice  of  a coordinate  system,  and  (2)  the  method  for  incrementing  the 
coordinates.  The  guiding  principle  in  making  these  two  choices  is  to 
try  to  equalize  the  contributions  to  the  integral  corresponding  to  each 
discrete  combination  of  the  variables  of  integration.  Unfortunately, 
different  source  locations  and  different  atmospheric  configurations 
require  somewhat  different  numerical  treatments  if  one  desires  to  mini- 
mize computational  expense.  This  led  to  the  development  of  two  separate 
computer  programs,  each  specialized  to  a particular  subset  of  the  possible 
configurations  of  the  problem.  The  differences  between  these  two  programs 
will  be  discussed  below. 

The  choice  of  an  appropriate  coordinate  system  for  the  integrals  is 
influenced  somewhat  by  the  method  chosen  for  incrementing  the  variables 
of  integration.  As  described  in  Section  7.1  equal  increments  were  chosen 
for  the  <J) ' integration.  Another  incrementing  scheme  which  is  relatively 
straightforward  to  implement  is  that  of  geometrically  increasing  (or 
decreasing)  increments.  According  to  this  scheme  the  coordinates  are 
incremented  in  a multiplicative,  rather  than  additive,  manner: 
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equal  ..ncrements:  z + Az  = z + Az 

geometric  increments:  z + Az'Otz  , a 1 

z z 

Each  new  discrete  value  of  z in  this  example  is  obtained  by  multiplying 

the  previous  discrete  value  by  a , hence  Az  = (a  -l)z  varies  in  direct 

z z 

proportion  to  the  value  of  z.  With  equal  increments  Az  would  be  constant. 
Of  course  in  the  z integral  required  to  calculate  E the  range  of  inte- 
gration goes  from  z = 0 to  z = 00 . With  either  incrementing  scheme  one 
would  need  an  infinite  number  of  increments  to  cover  this  range.  However, 
the  atmosphere  is  essentially  non-existent  above  50  km  altitude,  so  this 
can  be  used  as  an  upper  bo_nd  on  the  range  of  integration.  With  geo- 
metric increments  this  still  implies  that  an  infinite  number  of  increments 
are  needed.  The  technique  is  then  to  sum  the  contributions  to  the  inte- 
gral, starting  ct  the  upper  bound  (e.g.  40  km),  and  to  increment  the 
variable  of  integration  backwards: 

z - Az  = z/a  , a > 1 (92) 

z z 

The  running  total  of  these  contribute  ns  is  then  compared  to  each  new 
contribution  to  determine  when  this  total  has  more  or  less  converged. 
(Since  Az  keeps  getting  smaller  and  smaller,  eventually  a point  is 
reached  at  which  the  integrand,  multiplied  by  Az,  keeps  getting  smaller 
and  smaller  as  well  — except  for  extreme  cases  in  which  the  integral 
would  not  have  converged  analytically  either.  Such  is  never  the  case 
for  this  problem.)  When  convergence  is  detected,  an  analytic  approxima- 
tion to  che  remainder  of  the  integral  is  added  in  to  finish  the  compu- 
tation. In  this  case,  once  convergence  is  sensed,  the  integral  is  assumed 
to  be  constant  for  all  remaining  values  of  z. 

Since  the  atmosphere  is  ordinarily  densest  close  to  the  ground,  an 
important  portion  of  the  scattering  of  che  light  from  the  source  occurs 
at  low  altitudes.  Hence,  the  treatment  of  the  z integral  described  above 
was  found  to  be  both  useful  and  economical  for  some  configurations  of  the 
problem.  In  particular  one  of  the  two  computer  programs  written  to  solve 
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this  problem  uses  the  above  method  and  has  been  found  to  be  especially 
useful  when  the  source  is  close  to  the  ground,  outside  most  of  the 
earth's  atmosphere,  or  comfortably  above  a fog  layer.  However,  if  the 
source  is  close  to  a fog  layer,  immersed  in  a fog,  or  in  the  middle 
of  a hazy  atmosphere,  the  major  portion  of  the  singly  scattered  irradia- 
tion comes  from  scattering  centers  which  are  close  to  the  source  and 
which  consequently  receive  more  intense  radiation  from  it  than  scattering 
centers  farther  away.  In  these  latter  cases  one  needs  to  have  the 
smaller  increments  of  the  integration  occur  close  to  the  source.  In  this 
instance  the  z integration  is  separated  into  two  parts: 

z + Az  = z - (z  -z)/a  , 0 <.  t < z (93) 

s s z s 

2 - Az  = z + (z-z  )/a  , z < z <.  (z  +z  ) (94) 

s s z s s max7 

This  latter  z integration  scheme  is  used  in  the  other  of  the  two  com- 
puter programs  written  to  solve  this  problem.  The  method  used  in  this 
second  program  to  terminate  the  integral  series  is  the  same  as  that 

described  above.  It  has  been  found  that  values  of  a between  1.05  and 

z 

1.15,  corresponding  to  47.2  and  16.5  iterations  per  decade  (order  of 
magnitude),  respectively,  in  general  produce  results  for  the  z integral 
which  are  accurate  to  within  1%  to  5%.  This  accuracy  improves  for 
atmospheres  which  have  gentle  variations  in  density  as  a function  of 
altitude  (e.g.  hazes)  . The.  above  estimates  of  accuracy  are  based  on 
results  from  integration  by  the  trapezoidal  rule  (by  which  the  values 
of  the  integrand  at  the  endpoints  of  each  It.  interval  are  averaged  and 
are  then  multiplied  by  Az)  . ihe  trapezoidal  rule  is  used  for  all  the 
integrals  in  both  of  the  computer  programs  written  for  this  problem.  The 
different  methods  for  handling  the  integration  with  respect  to  z (alti- 
tude) are  the  only  difference  between  these  two  programs. 

Since  the  scattering  phase  functions  used  in  this  problem  tend  to 
scatter  major  amounts  of  radiation  in  the  forward  and  backward  directions, 
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much  of  the  singly  scattered  radiation  tends  to  follow  paths  close  to  a 
line  connecting  the  source  with  the  irradiated  spot  on  the  ground  which 
is  of  interest.  Consequently,  it  is  desirable  to  use  an  angular  coordi- 
nate system  for  the  remaining  integrals  which  facilitates  increasing  the 
number  of  increments  when  the  variables  of  Integration  are  close  to  the 
source  coordinates.  This  is  most  easily  accomplished  if  the  source  (or 
the  irradiated  spot)  is  near  an  origin  or  a ’’pole"  of  the  coordinate 
system.  Since  it  is  desirable  to  have  a handle  on  the  singly  scattered 
radiance  received  at  the  irradiated  spot  as  a function  of  the  angle  of 
incidence,  a spherical  coordinate  system  was  chosen  which  placed  the 
irradiated  spot  at  the  origin  and  placed  the  source  at  the  "pole".  This 
is  equivalent  to  rotating  a standard  spherical  coordinate  system  with  its 
"pole"  at  the  zenith  through  the  angle  a defined  in  Section  7.1.  The 
coordinates  in  this  rotated  system  are  designated  by  <p'  and  y'. 

Since  y'  is  the  cosine  of  an  angle  (6')  , equal  increments  in  y' 
correspond  to  unequal  increments  in  3'.  In  fact,  equal  increments  in  y' 
would  cause  the  smallest  increments  in  9'  to  occur  avound  6'  = 90°,  which 
is  just  the  opposite  of  what  is  needed.  The  next  logical  choice  might 
be  between  using  geometric  increments  in  6'  or  using  geometric  increments 
in  y'.  Geometric  increments  in  9'  would  prevent  the  use  of  the  trigo- 
nometric identities  used  for  the  4> 1 integral  (Section  7.1)  and  would 
require  that  sin0'  and  cos0'  be  recalculated  after  each  increment  in  0'. 
On  the  other  hand,  geometric  increments  in  u ' only  require  that  one 
square  root  be  calculated  after  each  increment: 

sin0 ' * \/ 1 - y ’ ^ (see  equation  for  y in  Section  7.1) 

This  latter  approach  is  used  in  loth  computer  programs,  in  which  case: 

y'  + Ay*  * 1 - (l-y’)/a^,  , ayl  > 1 
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A value  of  ay,  * 1.2,  corresponding  to  12.6  iterations  per  decade  (order 
of  magnitude),  has  been  found  to  produce  accuracies  within  10%  for  the  y' 
integral.  The  value  of  o , is  the  most  critical  for  the  contributions 
to  the  integral  near  y'  - 0 (01  = 90°),  since  the  first  effective  incre- 
ment in  0'  is  a.-,  large  as  A0'  ■ 9.6°  for  a . ■=  1.2. 

y 

When  tue  y*  integral  begins  to  converge,  the  remainder  of  the 
integral  is  computed  by  assuming  that  the  integrand  is  constant  for  the 
remaining  values  of  y',  and  this  remainder  is  added  to  the  running  sum 
to  obtain  the  total  integral.  One  needs  to  be  careful  with  this  approxi- 
mation. If  the  source  is  above  a low  lying  fog  layer,  unusually  large 
contributions  to  the  y'  integral  come  from  directions  which  are  somewhat 
distant  from  the  source  direction.  This  happens  because  there  are  a great 
many  scattering  centers  along  those  directions  which  cover  large  distances 
within  the  fog.  As  y*  approaches  unity  (01  approaches  0°)  the  distance 
for  which  these  rays  of  integration  remain  in  the  fog  layer  decreases 
dramatically.  At  the  same  time  the  scattering  phase  angle  falls  within 
the  range  x ^ 90°  (e.g.  45°  < x < 90°)  for  which  little  scattering  occurs. 
This  combination  of  effects  can  cause  the  y'  integral  to  appear  to  con- 
verge at  values  of  y'  corresponding  to  6'  ^ 50°.  Since  such  angles  are 
much  too  large  for  an  accurate  estimate  for  the  remainder  of  the  integral 
to  be  made,  one  must  input  to  the  computer  programs  an  unusually  stringent 
convergence  criterion  for  this  situation.  As  a check  on  this  problem,  the 
computer  programs  tabulate  the  z integral,  the  changes  to  the  y'  integral, 
and  the  running  sum  of  the  y 1 integral  for  each  value  of  y ' . 
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8 

CONCLUSIONS  ArTD  RECOMMENDATIONS 


In  this  report  we  have  developed  practical  atmospheric-radiative- 
transfer  models  which  are  used  to  ualculate  the  natural  and  artificial 
radiation  fields  in  atmospheres  ranging  in  visibility  from  the  clearest 
possible  -osphere  to  the  densest  fog.  The  model  for  calculation  of 
the  natural  illumination  due  to  solar  or  lunar  sources  is  valid  for  the 
spectral  range  M).30  pm  to  3.0  pm  in  that  multiple  scattering  by  gases 
and  aerosols  is  accounted  for.  At  the  present  time,  gaseous  absorption 
is  excluded  from  the  computational  model  but  absorption  by  particulates 
is  included.  It  is  also  assumed  that  the  surface  is  a uniform,  flat, 
Lambertian  (perfectly  diffuse)  reflector.  We  have  seen  that  the  spectral 
irradiance  (or  illuminance)  on  the  surface  is  dependent  upon  the  solar 
(or  lunar)  zenith  angle,  atmospheric  visibility,  end  surface  reflectance 
in  a somewhat  complicated  way.  The  surface  illumination  may  increase  or 
decrease  with  increasing  atmospheric  turbidity  depending  upon  the  surface 
reflectance.  These  effects  are  not  always  made  evident  when  using 
empirical  models  of  atmospheric  radiation  due  to  their  limited  range  of 
parameters.  Sound  physical  principles  have  been  used  in  the  development 
of  this  radiative-transfer  model  and  the  results  of  calculations  based 
upon  this  model  have  been  compared  with  results  from  more  elaborate  but 
also  more  lengthy  mathematical  calculations  and  the  agreement  is  excellent. 

A detailed,  mathematically  exact  model  has  also  been  developed  in 
this  report  which  was  used  to  calculate  the  singly-scattered  radiation 
field  in  any  uniform  atmosphere  due  to  point  sources  (flares).  Here,  we 
found  some  interesting  and  unexpected  results.  For  uniform  fogs,  a flare 
within  the  fog  produces  a highly  concentrated  distribution  of  radiation 
on  the  surface  whereas  a flare  above  the  fog  layer  produces  a much  more 
extended  distribution  of  radiation.  Hence,  if  it  is  desirable  to  create 
a relatively  broad  illumination  level  over  an  extended  area  then  flares 
should  be  ignited  well  above  a fog.  Some  caution  should  be  urged  in  the 
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interpretation  of  these  results  however,  since  multiple  scattering 
effects  have  not  been  included  in  the  analysis.  Inclusion  of  multiple 
scattering  would  probably  tend  to  broaden  the  spatial  distribution  due 
to  the  fact  that  scattering  always  tends  to  redistribute  the  radiation 
field.  As  indicated  in  Section  7,  there  were  many  computational  problems 
which  had  to  be  overcome  5n  the  detailed  numerical  integration  used  in 
the  point  source  model.  The  accuracy  of  the  results  were  tested  by 
taking  a point  source  (the  sun)  at  an  altitude  of  91  million  miles  and 
comparing  cne  irradiance  with  the  irradiance  as  calculated  according  to 
the  usual,  verified,  multiple-scattering  model. 

It  is  difficult  lo  compare  model  calculations  with  experimental 
results  since  the  experimentalists  rarely  define  or  estimate  critical 
parameters.  Many  times  a report  will  read  "partly  cloudy"  or  "low 
visibility",  which  is  quite  vague  and  does  not  quantify  the  necessary 
parameters  of  visual  range  or  optical  thickness.  What  is  needed  is  a 
detailed,  controlled  experiment  performed  under  a variety  of  environ- 
mental conditions.  More  reascnable  comparisons  can  then  be  made. 

Thus,  in  conclusion,  we  have  a multiple-scattering  model  which  can 
be  used  to  calculate  the  natural  illumination  level  in  clear,  hazy,  or 
foggy  atmospheres  as  a result  of  the  known  solar  or  lunar  input  irradiance 
at  the  top  of  the  atmosphere.  Also,  we  have  a single-scattering  point 
source  model  which  can  be  used  to  calculate  the  illumination  level  on 
the  surface  as  a result  of  any  number  of  flares  embedded  in  clear,  hazy, 
or  foggy  atmospheres . 

Before  a complete  model  of  the  radiation  field  can  be  made  for  all 
environmental  conditions,  there  is  some  additional  work  which  must  be 
done.  We  recommend  the  following: 

1.  Fog  and  Cloud  Model  Parameterization:  We  have  generated  the 

mathematical  methods  which  can  be  used  to  calculate  the  radiation  field 
in  an  atmosphere  in  which  flares  are  the  sources.  As  an  illustration  of 
the  computational  techniques  we  have  calculated  the  illumination  on  the 
surface  for  a simple  fog.  It  remains  to  be  seen  however,  to  what  extent 
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the  spatial  distribution  changes  if  different  fogs  or  cloud  layers  are 
present.  Perhaps  a flare  within  or  above  a cloud  layer  will  produce  a 
more  extended  spatial  distribution  of  illumination  than  one  within  a 
dense  fog.  Also,  multiple  scattering  should  be  included  in  the  point 
source  analysis,  especially  for  dense  fogs  or  clouds. 

2.  Temporal  Effects:  Specific  requirements  usually  require  that 

results  on  the  signatures  of  various  objects  in  a given  scene  be  generated 
within  minutes  or  hours  rather  than  days . For  this  reason , one  should 
consider  the  temporal  variation  in  the  environment.  For  example,  in 
thirty  minutes  a fog  can  change  its  optical  properties  considerably, 

and  meteorological  input  daU?  must  be  updated  in  a much  shorter  time 
interval.  Also,  wind-blown  dust  and  smoke  can  greatly  alter  the  illumi- 
nation level  for  a scene  in  a matter  of  seconds.  It  would  be  desirable, 
therefore,  to  model  the  radiation  field,  both  natural  and  artificial,  as 
a function  of  the  time-varying  environmental  parameters . 

3.  Distributed  Sources:  One  of  the  most  important  sources  of 

radiation,  particularly  at  night,  is  the  diffuse  sky  radiation  arising 
from  distant  terrain  sources  such  as  city  lights . Mountains  and  large 
hills  can  also  reflect  solar  or  lunar  radiation.  As  a result  of  advanced 
mathematical  techniques  partially  used  in  the  solution  of  the  point  source 
problem,  we  have  developed  the  capability  of  modeling  a distributed  source 
of  radiation  such  as  that  arising  from  a collection  of  city  lights, 
illuminated  clouds,  hillsides,  fires,  etc.  Clearly,  this  modeling  capa- 
bility shoula  be  implemented  for  target  signature  analysis  studies, 
especially  for  cloudless,  moonless  nights. 

4.  Surface  Effects:  In  almost  all  radiative-transfer  analysis  one 

assumes  a perfectly  diffuse  (Lambertian)  for  mathematical  simplicity.  In 
actual  practice,  however,  such  a surface  is  never  realized.  In  fact, 
strong  goniometric  properties  exist  for  many  natural  and  man-made  surfaces 
and  such  properties  can  affect  the  directly  reflected  radiation  by  a 
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considerable  amount.  Also,  it  is  recommended  that  reflection  of  flare 
radiation  by  the  surface  be  Included  in  the  continuing  analysis  of  the 
radiation  field.  Since  the  reflectance  of  vegetation  and  other  natu~al 
materials  is  quite  high  (^50-60%)  in  the  near  infrared  this  effect  could 
be  of  major  importance. 

5.  Target  Illumination:  Actual  targets  are  not  horizontal,  Hat 

objects  lying  on  horizontal  surfaces  but  are  extended  objects  charac- 
terized by  many  panels  oriented  at  all  directions  with  respect  to  the 
surface.  The  question  arises,  what  is  the  illumination  on  each  panel 
and  what  is  the  luminance  from  each  panel  as  a result  of  the  natural  and 
artificial  radiation  fields?  Some  empirical  models  have  been  developed 
which  attempt  to  answer  this  q lestion  for  daytime  natural  sources  such 
as  sunlight.  What  is  needed  however,  ic  a complete,  unified  model  which 
can  be  used  to  predict  the  illumination  on  specific  targets  for  flares 
and  natural  sources  of  illumination  for  various  environmental  conditions. 

6.  Thermal  Infrared:  It  is  usually  assumed  that  scattering  is 

only  of  importance  in  the  visible  part  of  the  spectrum.  While  that 
assumption  is  generally  true  for  clear  or  hazy  atmospheres  it  is  cer- 
ta?nly  not  true  for  clouds,  fogs,  or  dust  and  smoke  plumes.  The  reason 
is  that  the  particles  are  of  the  same  order  of  magnitude  as  the  wave- 
length in  the  infrared  (^5  ym  - 12  ym)  and  scattering  can  be  of  con- 
siderable importance.  Hence,  one  should  consider  scattering  and  emission 
of  radiation  by  particulates  for  the  thermal  region. 

7.  Other  Effects:  Finally,  we  recommend  that  additional  subjects 

be  investigated.  One  of  these  is  the  illumination  due  to  realistic 
anisotropic  flares.  It  is  assumed  that  a flare  is  represented  by  an 
isotropic  point  source,  whereas,  in  reality  it  has  an  intrinsic  radiation 
pattern  pproximately  that  of  a dipole  field.  This  effect  could  be  of 
importance  as  we  have  seen  in  the  current  analysis.  Another  effect  is 
the  radiation  field  arising  from  searchlight  illuminated  clouds.  This 
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is  a distributed  source  which  was  included  in  an  above  section.  Gaseous 
absorption  can  also  be  included  in  the  infrared.  We  have  models  which 
will  determine  the  radiance  and  transmittance  through  atmospheres  with 
various  profiles  of  humidity,  pressure,  and  temperature. 

It  is  believed  by  the  authors  that  if  the  above  recommendations 
were  carried  out  a much  more  comprehensive  model  would  result  which  could 
be  used  in  target  signature  analysis  investigations. 
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APPENDIX 


As  has  been  shown  throughout  this  report  the  illumination  due  to 
natural  and  artificial  sources  of  radiation  depends  upon  many  parameters 
and  the  results  cannot  be  represented  by  a single  graph  or  chart.  Never- 
theless, here  we  shall  present  a simple  description  of  the  overall  procedure 
whereby  a user  can  obtain  reasonably  accurate  values  for  the  illumination 
on  the  Earth's  surface  without  having  to  go  through  all  of  the  computational 
details  included  in  this  report. 

The  spectral  irradiance  due  to  the  moon  is  given  by  Figs.  28,  29,  and  30. 
Any  visibility  condition  from  a visual  range  of  336  km  down  to  a visual  range 
of  2 km  can  be  determined  by  interpolation  and  by  referring  to  Figs.  32  and 
33.  If  visual  range  is  not  available  to  the  user  in  field  applications  then 
an  appropriate  parameter  might  be  the  spectral  transmittance,  usually 
obtainable  by  a laser  transmissometer . In  this  case,  the  lunar  irradiance 
values  can  be  extended  to  very  hazy  or  foggy  conditions  as  illustrated  in 
Figs.  34,  35,  and  36.  Corresponding  values  for  any  other  wavelength  can  be 
determined  by  multiplying  these  curves  by  the  lunar  curve  of  Figure  21. 
Likewise,  any  lunar  phase  may  be  considered  by  multiplying  the  curves  by 
the  phase  function  of  Figure  19. 

A similar  procedure  can  be  given  for  the  description  of  the  illumination 
due  to  point  sources  of  radiation  such  as  flares.  If  we  consider  a single 
flare  with  an  output  spectral  power  of  one  watt  per  micrometer,  then  the 
irradiance  for  a flare  of  any  power  can  be  found  by  simply  multiplying  the 
power  by  the  curves  of  Figs.  39  through  48.  Most  of  the  curves  are  normalized 
to  a wavelength  of  0.55  um  but  any  wavelength  may  be  considered  by  inter- 
polation from  Figs,  45,  46,  47,  48,  and  by  considering  the  intrinsic  spectral 
power  output  of  a partio  1 ar  flare.  The  intrinsic  spectral  power  outputs  of 
two  flares  are  represented  n Figs.  49  and  50. 

Ire*,  i combination  of  flares  in  an  arbitrary  configuration  the  total 
Irradianre  on  the  surface  can  be  found  by  simply  addiug  the  results  for  each 
flare  separately.  Since  this  procedure  is  time  consuming  to  carry  out  for 


116 


each  configuration,  we  have  used  a computer  plotting  technique.  The  results 
for  a few  simple  configurations  are  presented  in  Figs.  52  through  56. 

The  above  is  a simple  description  of  a method  by  which  one  can  obtain 
illumination  levels  on  the  ground  for  moonlight  and  artificial  flares.  The 
basic  parameters  which  are  required  are:  visual  range  (km)  or  spectral 
transmittance  (dimensionless),  surface  albedo  (dimensionless),  date,  time 
of  day,  flare  configuration,  and  type  of  fog.  All  of  these  parameters  except 
the  last  are  relatively  easy  to  specify.  We  have  given  results  in  this 
report  based  upon  ‘'standard"  or  normal  fog  conditions.  If  there  are  unusual 
conditions  present  then  a more  detailed  treatment  is  necessary.  The  pro- 
cedures for  calculating  these  detailed  parameters  are  given  in  Sections  3 
and  4. 
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